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1995  SUMMER  RESEARCH  EXTENSION  PROGRAM  (SREP)  MANAGEMENT  REPORT 


1.0  BACKGROUND 


Under  the  provisions  of  Air  Force  Office  of  Scientific  Research  (AFOSR)  contract  F49620-90-C- 
0076,  September  1990,  Research  &  Development  Laboratories  (RDL),  an  8(a)  contractor  in 
Culver  City,  CA,  manages  AFOSR’ s  Summer  Research  Program.  This  report  is  issued  in  partial 
fulfillment  of  that  contract  (CLIN  0003 AC). 

The  Summer  Research  Extension  Program  (SREP)  is  one  of  four  programs  AFOSR  manages 
under  the  Summer  Research  Program.  The  Summer  Faculty  Research  Program  (SFRP)  and  the 
Graduate  Student  Research  Program  (GSRP)  place  college-level  research  associates  in  Air  Force 
research  laboratories  around  the  United  States  for  8  to  12  weeks  of  research  with  Air  Force 
scientists.  The  High  School  Apprenticeship  Program  (HSAP)  is  the  fourth  element  of  the  Summer 
Research  Program,  allowing  promising  mathematics  and  science  students  to  spend  two  months  of 
their  summer  vacations  working  at  Air  Force  laboratories  within  commuting  distance  from  their 
homes. 

SFRP  associates  and  exceptional  GSRP  associates  are  encouraged,  at  the  end  of  their  summer 
tours,  to  write  proposals  to  extend  their  summer  research  during  the  following  calendar  year  at 
their  home  institutions.  AFOSR  provides  funds  adequate  to  pay  for  SREP  subcontracts.  In 
addition,  AFOSR  has  traditionally  provided  further  funding,  when  available,  to  pay  for  additional 
SREP  proposals,  including  those  submitted  by  associates  from  Historically  Black  Colleges  and 
Universities  (HBCUs)  and  Minority  Institutions  (Mis).  Finally,  laboratories  may  transfer  internal 
funds  to  AFOSR  to  fund  additional  SREPs.  Ultimately  the  laboratories  inform  RDL  of  their 
SREP  choices,  RDL  gets  AFOSR  approval,  and  RDL  forwards  a  subcontract  to  the  institution 
where  the  SREP  associate  is  employed.  The  subcontract  (see  Appendix  1  for  a  sample)  cites  the 
SREP  associate  as  the  principal  investigator  and  requires  submission  of  a  report  at  the  end  of  the 
subcontract  period. 

Institutions  are  encouraged  to  share  costs  of  the  SREP  research,  and  many  do  so.  The  most 
common  cost-sharing  arrangement  is  reduction  in  the  overhead,  fringes,  or  administrative  charges 
institutions  would  normally  add  on  to  the  principal  investigator’s  or  research  associate’s  labor. 
Some  institutions  also  provide  other  support  (e.g.,  computer  run  time,  administrative  assistance, 
facilities  and  equipment  or  research  assistants)  at  reduced  or  no  cost. 

When  RDL  receives  the  signed  subcontract,  we  fund  the  effort  initially  by  providing  90  %  of  the 
subcontract  amount  to  the  institution  (normally  $18,000  for  a  $20,000  SREP).  When  we  receive 
the  end-of-research  report,  we  evaluate  it  administratively  and  send  a  copy  to  the  laboratory  for  a 
technical  evaluation.  When  the  laboratory  notifies  us  the  SREP  report  is  acceptable,  we  release 
the  remaining  funds  to  the  institution. 


Introduction  - 1 


2.0  THE  1995  SREP  PROGRAM 


SELECTION  DATA:  A  total  of  719  faculty  members  (SFRP  Associates)  and  286  graduate 
students  (GSRP  associates)  applied  to  participate  in  the  1994  Summer  Research  Program.  From 
these  applicants  185  SFRPs  and  121  GSRPs  were  selected.  The  education  level  of  those  selected 
was  as  follows: 


1994  SRP  Associates,  by  Degree 

SFRP 

GSRP 

PHD 

MS 

MS 

BS 

179 

6 

52 

69 

Of  the  participants  in  the  1994  Summer  Research  Program  90  percent  of  SFRPs  and  25  percent 
of  GSRPs  submitted  proposals  for  the  SREP.  Ninety  proposals  from  SFRPs  and  ten  from  GSRPs 
were  selected  for  funding,  which  equates  to  a  selection  rate  of  54  %  of  the  SFRP  proposals  and  of 
34%  for  GSRP  proposals. 


1995  SREP:  Proposals  Submitted  vs.  Proposals  Se 

lected 

Summer 

1994 

Participants 

Submitted 

SREP 

Proposals 

SREPs 

Funded 

SFRP 

185 

167 

90 

GSRP 

121 

29 

10 

TOTAL 

306 

196 

100 

The  funding  was  provided  as  follows: 


Contractual  slots  funded  by  AFOSR  75 

Laboratory  funded  14 

Additional  funding  from  AFOSR  11 

Total  100 
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Six  HBCU/MI  associates  from  the  1994  summer  program  submitted  SREP  proposals;  six  were 
selected  (none  were  lab-funded;  all  were  funded  by  additional  AFOSR  funds). 


!  Proposals  Submitted  and  Selected,  by  Laboratory 

Selected 

Armstrong  Laboratory 

41 

19 

Arnold  Engineering  Development  Center 

12 

4 

Frank  J.  Seiler  Research  Laboratory 

6 

3 

33 

19 

Rome  Laboratory 

31 

13 

Wilford  Hall  Medical  Center 

2 

1 

62 

37 

TOTAL 

Note:  Phillips  Laboratory  funded  3  SREPs;  Wright  Laboratory  funded  11;  and  AFOSR 
funded  11  beyond  its  contractual  75. 


The  306  1994  Summer  Research  Program  participants  represented  135  institutions. 


Institutions  Represented  on  the  1994  SRP  and  1995  SREP 

Number  of  schools 
represented  in  the 
Summer  92  Program 

Number  of  schools 
represented  in 
submitted  proposals 

Number  of  schools 
represented  in 
Funded  Proposals 

135 

118 

73 

Forty  schools  had  more  than  one  participant  submitting  proposals. 


Proposals  Submitted  Per  School 


Number  of  Proposals 


■  Submitted 
EJ  Selected 
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The  selection  rate  for  the  78  schools  submitting  1  proposal  (68%)  was  better  than  those 
submitting  2  proposals  (61%),  3  proposals  (50%),  4  proposals  (0%)  or  5+  proposals  (25%). 
The  4  schools  that  submitted  5+  proposals  accounted  for  30  (15%)  of  the  196  proposals 
submitted. 

Of  the  196  proposals  submitted,  159  offered  institution  cost  sharing.  Of  the  funded  proposals 
which  offered  cost  sharing,  the  minimum  cost  share  was  $1000.00,  the  maximum  was 
$68,000.00  with  an  average  cost  share  of  $12,016.00. 


Proposals  and  Institution  Cost  Sharing 

Proposals 

Submitted 

Proposals 

Funded 

159 

82 

Without  cost  sharing 

37 

18 

Total 

196 

100 

The  SREP  participants  were  residents  of  41  different  states.  Number  of  states  represented  at 
each  laboratory  were: 


States  Represented,  by  Proposals  Submitted/Selected  per  Laboratory 

Proposals 

Submitted 

Proposals 

Funded 

21 

13 

Arnold  Engineering  Development  Center 

5 

2 

Frank  J.  Seiler  Research  Laboratory 

5 

3 

16 

14 

14 

7 

Wilford  Hall  Medical  Center 

2 

1 

24 

20 

Eleven  of  the  1995  SEEP  Principal  Investigators  also  participated  in  the  1994  SREP. 


ADMINISTRATIVE  EVALUATION:  The  administrative  quality  of  the  SREP  associates’  final 
reports  was  satisfactory.  Most  complied  with  the  formatting  and  other  instructions  provided  to 
them  by  RDL.  Ninety  seven  final  reports  and  two  interim  reports  have  been  received  and  are 
included  in  this  report.  The  subcontracts  were  funded  by  $1,991,623.00  of  Air  Force  money. 
Institution  cost  sharing  totaled  $985,353.00. 
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TECHNICAL  EVALUATION:  The  form  used  for  the  technical  evaluation  is  provided  as 
Appendix  2.  ninety -two  evaluation  reports  were  received.  Participants  by  laboratory  versus 
evaluations  submitted  is  shown  below: 


Evaluations 

Percent 

Armstrong  Laboratory 

23 1 

20 

95.2 

Arnold  Engineering  Development  Center 

4 

4 

100 

Frank  J.  Seiler  Research  Laboratory 

3 

3 

100 

192 

18 

100 

Rome  Laboratory 

13 

13 

100 

Wilford  Hall  Medical  Center 

1 

1 

100 

37 

34 

91.9 

Total 

Notes: 

1:  Research  on  two  of  the  final  reports  was  incomplete  as  of  press  time  so  there  aren’t  any  technical 
evaluations  on  them  to  process,  yet.  Percent  complete  is  based  upon  20/21  =95.2% 

2:  One  technical  evaluation  was  not  completed  because  one  of  the  final  reports  was  incomplete  as  of 
press  time.  Percent  complete  is  based  upon  18/18=100% 

3:  See  notes  1  and  2  above.  Percent  complete  is  based  upon  93/97  =  95.9% 

The  number  of  evaluations  submitted  for  the  1995  SREP  (95.9%)  shows  a  marked 
improvement  over  the  1994  SREP  submittals  (65%). 

PROGRAM  EVALUATION:  Each  laboratory  focal  point  evaluated  ten  areas  (see  Appendix 
2)  with  a  rating  from  one  (lowest)  to  five  (highest).  The  distribution  of  ratings  was  as  follows: 


Rating 

Not  Rated 

1 

2 

3 

4 

5 

O^a! » If 

7 

1 

7 

kaumam 

The  8  low  ratings  (one  1  and  seven  2’s  )  were  for  question  5  (one  2)  “The  USAF  should 
continue  to  pursue  the  research  in  this  SREP  report”  and  question  10  (one  1  and  six  2’s)  “The 
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one-year  period  for  complete  SREP  research  is  about  right”,  in  addition  over  30%  of  the 
threes  (20  of  62)  were  for  question  ten.  The  average  rating  by  question  was: 


Question 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

4.6 

4.6 

4.7 

mm 

4.6 

4.7 

4.8 

4.5 

4.6 

4.0 

The  distribution  of  the  averages  was: 


AREA  AVERAGES 


hi  m 

IK 

II 1 

1 

1 

1 

1 

II  1 

1 

1 

1 

1 

4 

4.1 

4.2 

4.3 

4.4 

4.5 

4.6 

4.7 

4.8 

4.9 

5 

Area  10  “the  one-year  period  for  complete  SREP  research  is  about  right”  had  the  lowest 
average  rating  (4.1).  The  overall  average  across  all  factors  was  4.6  with  a  small  sample 
standard  deviation  of  0.2.  The  average  rating  for  area  10  (4.1)  is  approximately  three  sigma 
lower  than  the  overall  average  (4.6)  indicating  that  a  significant  number  of  the  evaluators  feel 
that  a  period  of  other  than  one  year  should  be  available  for  complete  SREP  research. 
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The  average  ratings  ranged  from  3.4  to  5.0.  The  overall  average  for  those  reports  that  were 
evaluated  was  4.6.  Since  the  distribution  of  the  ratings  is  not  a  normal  distribution  the  average 
of  4.6  is  misleading.  In  fact  over  half  of  the  reports  received  an  average  rating  of  4.8  or 
higher.  The  distribution  of  the  average  report  ratings  is  as  shown: 


It  is  clear  from  the  high  ratings  that  the  laboratories  place  a  high  value  on  AFOSR’s  Summer 
Research  Extension  Programs. 
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3.0  SUBCONTRACTS  SUMMARY 


Table  1  provides  a  summary  of  the  SREP  subcontracts.  The  individual  reports  are  published 
volumes  as  shown: 


Laboratory  Volume 

Armstrong  Laboratory  1A,  IB 

Arnold  Engineering  Development  Center  5 

Frank  J.  Seiler  Research  Laboratory  5 

Phillips  Laboratory  2 

Rome  Laboratory  3 

Wilford  Hall  Medical  Center  5 

Wright  Laboratory  4 A,  4B 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author  Author's  Sponsoring  n  .  Contract  Univ.  Cost 

Author’s  University _ Degree _ Lab  Performance  Period  Amount  Share 


Anderson ,  Janies  PhD 

Analytical  Chemistry  95-0807 

University  of  Georgia,  Athens,  GA 

Ashrafiuon ,  Hashem  PhD 

Mechanical  Engineering  95-0800 

Villanova  University,  Villanova,  PA 

Burke ,  Michael  PhD 

Tulane  University  95-0811 

Tulane  University,  New  Orleans,  LA 

Edwards ,  Paul  PhD 

Chemistry  95-0808 

Edinboro  Univ  of  Pennsylvania,  Edinboro,  PA 

Gerstman ,  Bernard  PhD 

Physics  95-0815 


Florida  International  Universi,  Miami,  FL 


Graetz ,  Kenneth  PhD 

Department  of  Psychology  95-0812 

University  of  Dayton,  Dayton,  OH 

Gupta ,  Pushpa  PhD 

Mathematics  95-0802 

University  of  Maine,  Orono,  ME 

Koch ,  Manfred  PhD 

Geophysics  95-0809 


Florida  State  University,  Tallahassee,  FL 


AL/EQ  01/01/95  12/31/95  $25000.00  $1826.00 

Determination  of  the  Redox  Capacity  of  Soil 
Sediment  and  Prediction  of  Pollutant 

AL/CF  01/01/95  12/31/95  $25000.00  $19528.00 

Finite  Element  Modeling  of  the  Human  Neck  and 
Its  Validation  for  the  ATB  Model 

AL/HR  01/01/95  09/30/95  $25000.00  $1818.00 

An  Examination  of  the  Validity  of  the  New  Air 
Force  ASVAB  Composites 

AL/EQ  01/01/95  12/31/95  $25000.00  $5000.00 

Fuel  Identification  by  Neural  Networks  Analysis 
of  the  Response  of  Vapor  Sensiti 

AL/OE  01/01/95  12/31/95  $24289.00  $2874.00 

A  Comparison  of  Multistep  vs  Singlestep 
Arrhenius  Integral  Models  for  Describing 

AL/HR  01/01/95  12/31/95  $25000.00  $0.00 

Effects  of  Mental  Workload  and  Electronic 
Support  on  Negotiation  Performance 

AL/AO  01/01/95  12/31/95  $25000.00  $2859.00 

Regression  to  the  Mean  in  Half  Life  Studies 


AL/EQ  12/01/94  04/30/95  $25000.00  $0.00 

Application  of  the  MT3D  Solute  Transport  Model 
to  the  Made-2  Site:  Calibration 


Novotny ,  Mark  PhD 

Supercomputer  Comp  Res.  I  95-0810 

Florida  State  University,  Tallahassee,  FL 


AL/EQ  01/01/95  12/31/95  $25000.00  $0.00 

Computer  Calculations  of  Gas-Phase  Reaction  Rate 

Constants 


AL/CF  01/01/95  12/31/95  $25000.00  $20550.00 

Surface  Fitting  Three  Dimensional  Human  Head 
Scan  Data 


Nurre ,  Joseph  PhD 

Mechanical  Engineering  95-0804 

Ohio  University,  Athens,  OH 

Piepmeier ,  Edward  PhD 

Pharmaceutics  95-0801 

University  of  South  Carolina,  Columbia,  SC 

Quinones ,  Miguel  PhD 

Psychology  95-0813 

Rice  University,  Houston,  TX 

Riccio ,  Gary  PhD 

Psychology  95-0806 

Univ  of  EL  Urbana-Champaign,  Urbana,  IL 

Shebilske ,  Wayne  PhD 

Dept  of  Psy  chology  95-0814 


Texas  A&M  University,  College  Station,  TX 


AL/AO  01/01/95  12/31/95  $25000.00  $11740.00 

The  Effects  of  Hyperbaric  Oxygenation  on 
Metabolism  of  Drugs  and  Other  Xenobioti 


AL/HR  01/01/95  12/31/95  $25000.00  $4000.00 

Maintaining  Skills  After  Training:  The  Role  of 
Opportunity  to  Perform  Trained  T 

AL/CF  01/01/95  05/31/95  $22931.00  $0.00 

Nonlinear  Transcutaneous  Electrical  Stimulation 
of  the  Vestibular  System 

AL/HR  01/01/95  12/31/95  $25000.00  $5614.00 

Cognitive  Factors  in  Distr  Training  Effects 

During  Acquisition  of  Complex  Skills 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author  Author's  Sponsoring  Contract  Univ.  Cost 

Author's  University _ Degree  Lab  Performance  Period  Amount  Share 

Weisenberger ,  Janet  PhD  AL/CF  01/01/95  12/31/95  S25000.00  S12234  00 

Dept  of  Speech*  Hearing  95-0805  Tactile  Feedback  for  Simulation  of  Object  Shape 

Ohio  State  University,  Columbus,  OH  and  Textural  Information  in  Hapt 

Hughes,  Rod  MA  AL/CF  01/01/95  12/31/95  S25000.00  S0.00 

Psychology  95-0803  Melatonin  Induced  Prophylactic  Sleep  as  a 

Oregon  Health  Sciences  University,  Portland,  OR  Countermeasure  for  Sleep  Deprivation 


Bapty ,  Theodore  MS  AEDC/E  01/01/95  12/31/95  S24979.00  SO  00 

Electrical  Engineering  95-0848  Plant-Wide  Preventive  Maintenance  &  Monitoring 

Vanderbilt  University,  Nashville,  TN 


Dorgan ,  John  PhD  FJSRL/F  01/01/95  12/31/95  $25000.00  $0.00 

Chemical  Engineering  95-0834  Block  Copolymers  at  Inorganic  Solid  Surfaces 

Colorado  School  of  Mines,  Golden,  CO 


Jungbauer ,  Mary  Ann  PhD 

Chemistry  95-0836 

Barry  University,  Miami,  FL 

Statman  ,  David  PhD 

Physics  95-0835 

Allegheny  College,  Meadville,  PA 

,  Krishnaswamy  PhD 

Aeronautics  95-0818 

University  of  Houston,  Houston,  TX 

Ashgriz ,  Nasser  PhD 

Mechanical  Engineering  95-0816 

SUNY-Buffalo,  Buffalo,  NY 

Bellem ,  Raymond  PhD 

Computer  Science  95-0817 

Embry-Riddle  Aeronautical  Univ,  Prescott,  A Z 

Brzosko ,  Jan  PhD 

Nuclear  Physics  95-082  8 

Stevens  Institute  of  Tech,  Hoboken,  NJ 

Damodaran ,  Meledath  PhD 

Math  &  Computer  Science  95-0831 


University  of  Houston-Victoria,  Victoria,  TX 


DeLyser ,  Ronald  PhD 

Electrical  Engineering  95-0877 

University  of  Denver,  Denver,  CO 


Diels ,  Jean-Claude  PhD 

Physics  95-0819 

University  of  New  Mexico,  Albuquerque,  NM 


Henson  ,  James  PhD 

Electrical  Engineering  95-082  0 

University  of  Nevada,  Reno,  NV 


Kaiser ,  Gerald  PhD 

Physics  95-0821 

University  of  Mass/Lowell,  Lowell,  MA 


FJSRL/F  01/01/95  12/31/95  $25000.00  $24714.00 

Non-Linear  Optical  Properties  of  Polyacetylenes 

and  Related  Substituted  Compound 

FJSRL/F  01/01/95  12/31/95  $25000.00  $6500.00 

Studies  of  Second  Harmonic  Generation  in  Glass 
Waveguides 

PL/RK  01/01/95  12/31/95  $24993.00  $8969.00 

Mixed-Mode  Fracture  of  Solid  Propellants 


PL/RK  01/01/95  12/31/95  $25000.00  $22329.00 

Effects  of  the  Jet  Characteristics  on  the 
Atomization  and  Mixing  in  A  Pair  of  Im 

PL/VT  12/01/94  11/30/95  $20000.00  $8293.00 

Experimental  Studies  of  the  Effects  of  Ionizing 
Radiation  on  Commerically  Proces 

PL/WS  11/01/94  02/01/95  $24943.00  $0.00 

Neutron  Diagnostics  for  Pulsed  Plasmas  of 
Compact  Toroid  -  Marauder  Type 

PL/LI  01/01/95  12/31/94  $24989.00  $9850.00 

Parallel  Computation  of  Zernike  Aberration 
Coefficients  for  Optical  Aber  Correct 

PL/WS  01/01/95  12/31/95  $25000.00  $46066.00 

Quality  Factor  Evaluation  of  Complex  Cavities 

PL/LI  01/01/95  12/31/95  $25000.00  $0.00 

Unidirectional  Ring  Lasers  and  Laseer  Gyros  with 

Multiple  Quantum  Well  Gain  Medi 

PL/WS  01/01/95  12/31/95  $25000.00  $0.00 

Automatic  Feature  Extraction  and  Assessment  of 
Wideband  Range-Doppler  Imagery  of 

PL/GP  01/01/95  12/31/95  $25000.00  $5041.00 

Multiresolution  Analysis  with  Physical  Wavelets 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author  Author's 

Author's  University _ Degree 

KowaJak ,  Albert  PhD 

Chemistry  95-0822 

University  of  Massachusetts/Lo,  Lowell,  MA 


Malloy ,  Kevin  PhD 

Electrical  Engineering  95-0829 

University  of  New  Mexico,  Albuquerque,  NM 

Prasad ,  Sudhakar  PhD 

Physics  95-0823 

University  of  New  Mexico,  Albuquerque,  NM 

Purtill ,  Mark  PhD 

Mathematics  95-0824 

Texas  A&M  Univ-Kingsville,  Kingsville,  TX 

Rudolph ,  Wolfgang  PhD 

Physics  95-0833 

University  of  New  Mexico,  Albuquerque,  NM 

Stone ,  Alexander  PhD 

Mathematics  &  Statistics  95-0827 

University  of  New  Mexico,  Alburquerque,  NM 

Swenson ,  Charles  PhD 

Dept  of  Electrical  Engnr  95-0826 

Utah  State  University,  Logan,  UT 

Lipp ,  John  MS 

Electrical  Engineering  95-0832 

Michigan  Technological  Univ,  Houghton,  MI 

Petroski ,  Janet  BA 

Chemistry  95-0830 

Cal  State  Univ/Northridge,  Northridge,  CA 

Salasovich ,  Richard  MS 

Mechanical  Engineering  95-0825 

University  of  Cincinnati,  Cincinnati,  OH 

Aalo ,  Valentine  PhD 

Dept  of  Electrical  Engr  95-0837 

Florida  Atlantic  University,  Boca  Raton,  FL 

Amin ,  Moeness  PhD 

Electrical  Engineering  95-0838 

Villanova  University,  Villanova,  PA 

Benjamin ,  David  PhD 

Computer  Science  95-0839 

Oklahoma  State  University,  Stillwater,  OK 

Choudhury ,  Ajit  PhD 

Engineering  95-0840 

Howard  University,  Washington,  DC 

Harackiewicz ,  Frances  PhD 

Electrical  Engineering  95-0841 

So.  Dlinois  Univ-Carbondale,  Carbondale,  IL 


Sponsoring  „  ,  _  .  .  Contract  Univ.  Cost 

Lab  Performance  Penod  Ammnt  Share 


PL/GP  01/01/95  12/31/95  S24996.00  S4038.00 

The  Synthesis  and  Chemistry  of  Peroxonitrites 
and  Peroxonitrous  Acid 

PL/VT  01/01/95  12/31/95  S24999.00  S0.00 

Temperature  &  Pressure  Dependence  of  the  Band 
Gaps  &  Band  Offsets 

PL/LI  01/01/95  12/31/95  S25000.00  SI  1047.00 

Theoretical  Studies  of  the  Performance  of  Novel 
Fiber-Coupled  Imaging  Interferom 

PL/WS  01/01/95  12/31/95  S25000.00  S100.00 

Static  and  Dynamic  Graph  Embedding  for  Parallel 

Programming 

PL/LI  01/01/95  12/31/95  S24982.00  S6000.00 

Ultrafast  Process  and  Modulation  in  Iodine 

Lasers 

PL/WS  01/01/95  12/31/95  S24969.00  S0.00 

Impedance  Matching  And  Reflection  Minimization 
For  Transient  EM  Pulses  Through  D 

PL/VT  01/01/95  12/31/95  $25000.00  S25000.00 

Low  Power  Retromodulator  based  Optical 
Transceiver  for  Satellite  Communications 

PL/LI  01/01/95  12/31/95  S24340.00  S15200.00 

Improved  Methods  of  Tilt  Measurement  for 
Extended  Images  in  the  Presence  of  Atmo 

PL/RK  10/01/94  12/31/94  S4279.00  S0.0O 

Thermoluminescence  of  Simple  Species  in 
Molecular  Hydrogen  Matrices 

PL/VT  01/01/95  12/31/95  $25000.00  $4094.00 

Design,  Fabrication,  Intelligent  Cure,  Testing, 
and  Flight  Qualification  of  an  A 

RL/C3  01/01/95  12/31/95  $25000.00  $13120.00 

Performance  Study  of  an  ATM/Satellite  Network 

RL/C3  01/01/95  12/31/95  $25000.00  $34000.00 

Interference  Excision  in  Spread  Spectrum 
Communication  Systems  Using  Time-Freque 

RL/C3  01/01/95  12/31/95  $24970.00  $0.00 

Designing  Software  by  Decomposition  using  KIDS 

RL/OC  11/30/94  10/31/95  $25000.00  $0.00 

Detection  Performance  of  Over  Resolved  Targets 
with  Non-Uniform  and  Non-GAussian 

RL/ER  01/01/95  12/31/95  $23750.00  $29372.00 

Computer-Aided-Design  Program  for  Solderless 
Coupling  Between  Microstrip  and  Str 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author  Author's 

Author's  University _ Degree 

Losiewicz ,  Beth  PhD 

Psycholinguistics  95-0842 

Colorado  State  University,  Fort  Collins,  CO 

Musavi ,  Mohamad  PhD 

95-0843 

University  of  Maine,  Orono,  ME 

Norgard  ,  John  PhD 

Elec  &  Comp  Engineering  95-0844 

Univ  of  Colorado-Colorado  Sprg,  Colorado 

Richardson ,  Dean  PhD 

Photonics  95-0845 

SUNY  Institute  of  Technology,  Utica,  NY 

Ryder,  Jr. ,  Daniel  PhD 

Chemical  Engineering  95-0846 

Tufts  University,  Medford,  MA 


Sponsoring 

Lab 


,  Contract 

Performance  Period  Amount 


Univ.  Cost 
Share 


RL/IR  01/01/95  12/31/95  S25000.00  $4850.00 

Spanish  Dialect  Identification  Project 


RL/IR  01/01/95  12/31/95  S25000.00  S12473.00 

Automatic  Image  Registration  Using  Digital 
Terrain  Elevation  Data 

RL/ER  01/01/95  12/31/95  $25000.00  $2500.00 

Infrared  Images  of  Electromagnetic  Fields 


RL/OC  01/01/95  12/31/95  S25000.00  $15000.00 

Femtosecond  Pump-Probe  Spectroscopy  System 


RL/ER  01/01/95  12/31/95  S25000.00  $0.00 

Synthesis  and  Properties  B-Diketonate-Modified 
Heterobimetallic  Alkoxides 


Zhang ,  Xi-Cheng  PhD 

Physics  95-0847 

Rensselaer  Polytechnic  Institu,  Troy,  NY 

Drost-Hansen  ,  Walter  PhD 

Chemistry  95-0875 

University  of  Miami,  Coral  Gables,  FL 

Baginski ,  Michael  PhD 

Electrical  Engineering  95-0869 

Auburn  University,  Auburn,  AL 

Berdichevsky ,  Victor  PhD 

Aerospace  Engineering  95-0849 

Wayne  State  University,  Detroit,  MI 

Buckner ,  Steven  PhD 

Chemistry  95-0850 

Colullmbus  College,  Columbus,  GA 

Carroll ,  James  PhD 

Electrical  Engineering  95-0881 

Clarkson  University,  Potsdam,  NY 

Choate ,  David  PhD 

Mathematics  95-0851 

Transylvania  University,  Lexington,  KY 

Clarson  ,  Stephen  PhD 

Materials  Sci  &  Eng  95-0852 


University  of  Cincinnati,  Cincinnati,  OH 


Cone ,  Milton  PhD 

Comp  Science  &  Elec  Eng  95-0853 

Embry-Riddel  Aeronautical  Univ,  Prescott,  A Z 

Courier ,  Robert  PhD 

Mechanical  Engineering  95-0854 

Louisiana  State  University,  Baton  Rouge,  LA 


RL/ER  01/01/95  12/31/95  $25000.00  $0.00 

Optoelectronic  Study  of  Seniconductor  Surfaces 
and  Interfaces 

WHMC/  01/01/95  12/31/95  $25000.00  $8525.00 

Biochemical  &  Cell  Physiological  Aspects  of 
Hyperthermia 

WL/MN  01/01/95  12/31/95  $24995.00  S  10098.00 

An  Investigation  of  the  Heating  and  Temperature 
Distribution  in  Electrically  Exc 

WL/FI  01/01/95  12/31/95  $25000.00  $0.00 

Micromechanics  of  Creep  in  Metals  and  Ceramics 
at  High  Temperature 

WL/PO  01/01/95  12/31/95  $24900.00  $8500.00 

Development  of  a  Fluorescenece-Based  Chemical 
Sensor  for  Simultaneous  Oxygen  Qua 

WL/PO  01/01/95  12/31/95  $24944.00  $38964.00 

Development  of  High-Performance  Active 
Dynamometer  System  for  Machines  and  Drive 

WL/AA  01/01/95  12/31/95  $24993.00  $8637.00 

SOLVING  z (t)  =  ln{/A[cos (wit) ] +B[sin (w2t) ] +C) 

WL/ML  12/01/94  11/30/95  S25000.00  $15000.00 

Synthesis,  Processing  and  Characterization  of 
Nonlinear  Optical  Polymer  Thin  Fil 

WL/AA  01/01/95  12/31/95  $25000.00  $11247.00 

An  Investigation  of  Planning  and  Scheduling 
Algorithms  for  Sensor  Management 

WL/MN  01/01/95  12/31/95  $25000.00  $3729.00 

A  Study  to  Determine  Wave  Gun  Firing  Cycles  for 

High  Performance  Model  Launches, 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author  Author's 

Author's  University  _  Degree 

Dominic ,  Vincent  PhD 

Electro  Optics  Program  95-0868 

University  of  Dayton,  Dayton,  OH 

Fadel ,  Georges  PhD 

Dept  of  Mechanical  Engnr  95-0855 

Clemson  University,  Clemson,  SC 

Gould ,  Richard  PhD 

Mechanical  Engineering  95-0856 

North  Carolina  State  Univ,  Raleigh,  NC 


Sponsoring  .  _  .  Contract  Univ.  Cost 

Lab  Performance  Period  Amount  Share 

WL/ML  01/01/95  12/31/95  $25000.00  S12029.00 

Characterization  of  Electro-Optic  Polymers 


Hardie ,  Russell  PhD 

Electrical  Engineering  95-0882 

Univsity  of  Dayton,  Dayton,  OH 

Hodel ,  Alan  PhD 

Electrical  Engineering  95-0870 

Auburn  University,  Auburn,  AL 

Janus ,  Jonathan  PhD 

Aerospace  Engineering  95-0871 

Mississippi  State  University,  Mississippi  State, 

Jasiuk ,  Iwona  PhD 

Dept  of  Materials  Science  95-0857 

Michigan  State  University,  East  Lansing,  MI 

Jouny ,  Ismail  PhD 

Electrical  Engineering  95-0880 

Lafayette  College,  Easton,  PA 

Li ,  Jian  PhD 

Electrical  Engineering  95-0859 

University  of  Florida,  Gainesville,  FL 

Lin ,  Chun-Shin  PhD 

Electrical  Engineering  95-0883 

University  of  Missouri-Columbi,  Columbia,  MO 

Lin ,  Paul  PhD 

Mechanical  Engineering  95-0860 

Cleveland  State  University,  Cleveland,  OH 

Liou  ,  Juin  PhD 

Electrical  Engineering  95-0876 

University  of  Central  Florida,  Orlando,  FL 

Nandhakumar ,  Nagaraj  PhD 

Electrical  Engineering  95-0861 

University  of  Virginia,  Charlottesville,  VA 

Pasala ,  Krishna  PhD 

Dept  of  Electrical  Engr  95-0879 

University  of  Dayton,  Dayton,  OH 

Perkowski ,  Marek  PhD 

Dept  of  Electrical  Engnr  95-0878 

Portland  State  University,  Portland,  OR 


WL/MT  01/01/95  12/31/95  S25000.00  $8645.00 

A  Methodology  for  Affordability  in  the  Design 

Process 

WL/PO  01/01/95  12/31/95  $24998.00  $9783.00 

Data  Reduction  and  Analysis  for  laser  Doppler 
Velocimetry 

WL/AA  01/01/95  12/31/95  $24999.00  $7415.00 

Hyperspectral  Target  Identification  Using  Bomen 
Spectometer  Data 

WL/MN  01/01/95  12/31/95  $24990.00  $9291.00 

Robust  Falut  Tolerant  Control:  Fault  Detection 
and  Classification 

WL/MN  01/01/95  12/31/95  $25000.00  $7143.00 

Multidimensional  Algorithm  Development  & 
Analysis 

WL/ML  01/01/95  12/31/95  $25000.00  $0.00 

Characterization  of  Interfaces  in  Metal-Matrix 

Composites 

WL/AA  01/01/95  12/31/95  $24300.00  $5200.00 

TSI  Mitigation:  A  Mountaintop  Database  Study 

WL/AA  10/10/95  12/31/95  $25000.00  $4000.00 

Comparative  Study  and  Performance  Analysis  of 
High  Resolution  SAR  Imaging  Techni 

WL/FI  01/01/95  12/31/95  $25000.00  $2057.00 

Prediction  of  Missile  Trajectory 

WL/FI  01/01/95  12/31/95  $25000.00  $6886.00 

Three  Dimensional  Deformation  Comparison  Between 
Bias  and  Radial  Aircraft  Tires 

WL/EL  01/01/95  12/31/95  $25000.00  $11040.00 

Investigation  of  AlGaAs/GaAs  Heterojunction 
Bipolar  Transister  Reliability  Based 

WL/AA  01/01/95  12/31/95  $24979.00  $4500.00 

Thermophysical  Invariants  fro,  LWIR  Imagery  for 

ATR 

WL/AA  01/01/95  12/31/95  S25000.00  $1078.00 

Effect  of  Electromagmetic  Enviornment  on  Array 
Signal  Processing 


WL/AA  01/01/95  09/15/95  $24947.00  $1 

Functional  Decomposition  of 
Binary, Multiple-Valued,  &  Fuzzy  Logic 


$18319.00 
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1995  SREP  SUB-CONTRACT  DATA 


Report  Author 

Author's  University’ _ 

Reeves ,  Stanley 

Dept  of  Electrical  Engnr 

Auburn  University,  Auburn,  AL 


Author's 

Degree 


PhD 

95-0862 


Sponsoring  performance  perj0(j 


Contract 

Amount 


Univ.  Cost 
Share 


WL/MN  01/01/95  12/31/95  S25000.00  S0.00 

Superresolution  of  Passive  Millimeter-Wave 
Imaging 


Rule ,  William  PhD 

Engineering  Mechanics  95-0872 

University  of  Alabama,  Tuscaloosa,  AL 

Schauer ,  John  PhD 

Mech  &  Aerosp  Eng  95-0873 

University  of  Dayton,  Dayton,  OH 


WL/MN  01/01/95  12/31/95  S24968.00  $14576.00 

Development  of  a  Penetrator  Optimizer 


WL/PO  11/01/94  11/30/95  $25000.00  $7428.00 

Heat  Transfer  for  Turbine  Blade  Film  Cooling 
with  Free  Stream  Turbulence  -  Measu 


Schwartz ,  Carla  PhD 

Electrical  Engineering  95-0863 

University  of  Florida,  Gainesville,  FL 


WL/FI  01/01/95  12/31/95  $25000.00  $0.00 

Neural  Network  Identification  and  Control  in 
Metal  Forging 


Simon ,  Terrence  PhD 

Dept  of  Mechanical  Engineering  95-0864 
University  of  Minnesota,  Minneapolis,  MN 


WL/PO  01/01/95  12/31/95  $24966.00  $3996.00 

Documentation  of  Separating  and  Separated 
Boundary  Layer  Flow,  for  Application 


Skowronski ,  Marek  PhD 

Solid  State  Physics  95-0865 

Carnegie  Melon  University,  Pittsburgh,  PA 


WL/EL  01/01/95  12/31/95  $25000.00  $6829.00 

Transmission  Electron  Microscopy  of 
Semiconductor  Heterojunctions 


Thirunarayan  ,  Krishnaprasad  PhD 

Computer  Science  95-0866 

Wright  State  University,  Dayton,  OH 


WL/EL  01/01/95  12/31/95  $25000.00  $2816.00 

VHDL-93  Parser  in  SWI-PROLOG:  A  Basis  for  Design 

Query  System 


Trelease ,  Robert  PhD 

Dept  of  Anatomy  &  Cell  Bi  95-0867 

University  of  California,  Los  Angeles,  CA 


WL/ML  12/01/94  12/01/95  $25000.00  $0.00 

Development  of  Qualitative  Process  Control 
Discovery  Systems  for  Polymar  Composi 


Tsai ,  Chi-Tay  PhD 

Engineering  Mechanics  95-0874 

Florida  Atlantic  University  ,  Boca  Raton,  FL 


WL/MN  01/01/95  12/31/95  $24980.00  $0.00 

Improved  Algorithm  Development  of  Massively 
Parallel  Epic  Hydrocode  in  Cray  T3D 


Lewis ,  John  MS 

Materials  Science  Engrng  95-0858 

University  of  Kentucky,  Lexington,  KY 


WL/MN  01/01/95  12/31/95  $25000.00  $13833.00 

The  Characterization  of  the  Mechanical 
Properties  of  Materials  in  a  Biaxial  Stre 
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APPENDIX  1: 

SAMPLE  SREP  SUBCONTRACT 
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AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH 
1995  SUMMER  RESEARCH  EXTENSION  PROGRAM 
SUBCONTRACT  95-0837 

BETWEEN 


Research  &  Development  Laboratories 
5800  Uplander  Way 
Culver  City,  CA  90230-6608 

AND 

Florida  Atlantic  University 
Department  of  Electrical  Engineering 
Boca  Raton,  FL  3343 1 


REFERENCE:  Summer  Research  Extension  Program  Proposal  95-0837 

Start  Date:  01-01-95  End  Date:  12-31-95 

Proposal  Amount:  $25,000.00 

(I)  PRINCIPAL  INVESTIGATOR:  Dr  Valentine  A  Aalo 

Department  of  Electrical  Engineering 
Florida  Atlantic  University 
Boca  Raton,  FL  33431 


(2)  UNITED  STATES  AFOSR  CONTRACT  NUMBER:  F49620-93-C-0063 


(3)  CATALOG  OF  FEDERAL  DOMESTIC  ASSISTANCE  NUMBER  (CFD A):  12.800 
PROJECT  TITLE:  AIR  FORCE  DEFENSE  RESEARCH  SOURCES  PROGRAM 


(4) 


ATTACHMENT  1 
2 
3 


REPORT  OF  INVENTIONS  AND  SUBCONTRAACT 

CONTRACT  CLAUSES 

FINAL  REPORT  INSTRUCTIONS 


***SIGN  SREP  SUBCONTRACT  AND  RETURN  TO  RDL*** 
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1.  BACKGROUND:  Research  &  Development  Laboratories  (RDL)  is  under  contract 
(F49620-93-C-0063)  to  the  United  States  Air  Force  to  administer  the  Summer 
Research  Program  (SRP),  sponsored  by  the  Air  Force  Office  of  Scientific  Research 
(AFOSR),  Bolling  Air  Force  Base,  D.C.  Under  the  SRP,  a  selected  number  of  college 
faculty  members  and  graduate  students  spend  part  of  the  summer  conducting  research 
in  Air  Force  laboratories.  After  completion  of  the  summer  tour  participants  may 
submit,  through  their  home  institutions,  proposals  for  follow-on  research.  The  follow- 
on  research  is  known  as  the  Summer  Research  Extension  Program  (SREP). 
Approximately  61  SREP  proposals  annually  will  be  selected  by  the  Air  Force  for 
funding  of  up  to  $25,000;  shared  funding  by  the  academic  institution  is  encouraged. 
SREP  efforts  selected  for  funding  are  administered  by  RDL  through  subcontracts  with 
the  institutions.  This  subcontract  represents  an  agreement  between  RDL  and  the 
institution  herein  designated  in  Section  5  below. 

2.  RDL  PAYMENTS:  RDL  will  provide  the  following  payments  to  SREP  institutions: 

•  80  percent  of  the  negotiated  SREP  dollar  amount  at  the  start  of  the  SREP 
research  period. 

•  The  remainder  of  the  funds  within  30  days  after  receipt  at  RDL  of  the 
acceptable  written  final  report  for  the  SREP  research. 

3.  INSTITUTION’S  RESPONSIBILITIES:  As  a  subcontractor  to  RDL,  the  institution 
designated  on  the  title  page  will: 
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a.  Assure  that  the  research  performed  and  the  resources  utilized  adhere  to  those 
defined  in  the  SREP  proposal. 

b.  Provide  the  level  and  amounts  of  institutional  support  specified  in  the  SREP 
proposal.. 

c.  Notify  RDL  as  soon  as  possible,  but  not  later  than  30  days,  of  any  changes  in 
3a  or  3b  above,  or  any  change  to  the  assignment  or  amount  of  participation  of 
the  Principal  Investigator  designated  on  the  title  page. 

d.  Assure  that  the  research  is  completed  and  the  final  report  is  delivered  to  RDL 
not  later  than  twelve  months  from  the  effective  date  of  this  subcontract,  but  no 
later  than  December  31,  1998.  The  effective  date  of  the  subcontract  is  one 
week  after  the  date  that  the  institution’s  contracting  representative  signs  this 
subcontract,  but  no  later  than  January  15,  1998. 

e.  Assure  that  the  final  report  is  submitted  in  accordance  with  Attachment  3. 

f.  Agree  that  any  release  of  information  relating  to  this  subcontract  (news 
releases,  articles,  manuscripts,  brochures,  advertisements,  still  and  motion 
pictures,  speeches,  trade  associations  meetings,  symposia,  etc.)  will  include  a 
statement  that  the  project  or  effort  depicted  was  or  is  sponsored  by:  Air  Force 
Office  of  Scientific  Research,  Bolling  AFB,  D.C. 

g.  Notify  RDL  of  inventions  or  patents  claimed  as  the  result  of  this  research  as 
specified  in  Attachment  1 . 

h.  RDL  is  required  by  the  prime  contract  to  flow  down  patent  rights  and  technical 
data  requirements  to  this  subcontract.  Attachment  2  to  this  subcontract 
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contains  a  list  of  contract  clauses  incorporated  by  reference  in  the  prime 
contract. 


4.  All  notices  to  RDL  shall  be  addressed  to: 

RDL  AFOSR  Program  Office 
5800  Uplander  Way 
Culver  City,  CA  90230-6609 

5.  By  their  signatures  below,  the  parties  agree  to  provisions  of  this  subcontract. 


Abe  Sopher  Signature  of  Institution  Contracting  Official 

RDL  Contracts  Manager 


Typed/Printed  Name 


Date 


Title 


Institution 


Date/Phone 
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ATTACHMENT  2 
CONTRACT  CLAUSES 


This  contract  incorporates  by  reference  the  following  clauses  of  the  Federal  Acquisition 
Regulations  (FAR),  with  the  same  force  and  effect  as  if  they  were  given  in  full  text.  Upon 
request,  the  Contracting  Officer  or  RDL  will  make  their  full  text  available  (FAR  52.252-2). 

FAR  CLAUSES  TITLE  AND  DATE 


52.202- 1  DEFINITIONS 

52.203- 3  GRATUITIES 


52.203- 5 

52.203- 6 

52.203- 7 

52.203- 8 

52.203- 10 

52.203- 12 

52.204- 2 
52.209-6 


52.212-8 

52.215- 2 

52.215- 10 


COVENANT  AGAINST  CONTINGENT  FEES 

RESTRICTIONS  ON  SUBCONTRACTOR 
SALES  TO  THE  GOVERNMENT 

ANTI-KICKBACK  PROCEDURES 

CANCELLATION,  RECISSION,  AND 

RECOVERY  OF  FUNDS  FOR  ILLEGAL  OR  IMPROPER 

ACTIVITY 

PRICE  OR  FEE  ADJUSTMENT  FOR  ILLEGAL 
OR  IMPROPER  ACTIVITY 

LIMITATION  ON  PAYMENTS  TO  INFLUENCE 
CERTAIN  FEDERAL  TRANSACTIONS 

SECURITY  REQUIREMENTS 

PROTECTING  THE  GOVERNMENT’S 
INTEREST  WHEN  SUBCONTRACTING  WITH 
CONTRACTORS  DEBARRED,  SUSPENDED,  OR 
PROPOSED  FOR  DEBARMENT 

DEFENSE  PRIORITY  AND  ALLOCATION 
REQUIREMENTS 

AUDIT  AND  RECORDS  -  NEGOTIATION 

PRICE  REDUCTION  FOR  DEFECTIVE  COST 
OR  PRICING  DATA 
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SUBCONTRACTOR  COST  OR  PRICING 
DATA 


52.215- 14 

52.215- 8 
52.215.18 

52.222- 3 

52.222- 26 

52.222- 35 

52.222- 36 

52.222- 37 

52.223- 2 

52.223- 6 

52.224- 1 

52.224- 2 

52.225- 13 

52.227- 1 

52.227- 2 


INTEGRITY  OF  UNIT  PRICES 

ORDER  OF  PRECEDENCE 

REVERSION  OR  ADJUSTMENT  OF  PLANS 
FOR  POSTRETIREMENT  BENEFITS  OTHER 
THAN  PENSIONS 

CONVICT  LABOR 

EQUAL  OPPORTUNITY 

AFFIRMATIVE  ACTION  FOR  SPECIAL 
DISABLED  AND  VIETNAM  ERA 
VETERANS 

AFFIRMATIVE  ACTION  FOR 
HANDICAPPED  WORKERS 

EMPLOYMENT  REPORTS  ON  SPECIAL 
DISABLED  VETERAN  AND  VETERANS  OF  THE 
VIETNAM  ERA 

CLEAN  AIR  AND  WATER 

DRUG-FREE  WORKPLACE 

PRIVACY  ACT  NOTIFICATION 

PRIVACY  ACT 

RESTRICTIONS  ON  CONTRACTING  WITH 
SANCTIONED  PERSONS 

ALT.  I  -  AUTHORIZATION  AND  CONSENT 

NOTICE  AND  ASSISTANCE  REGARDING 
PATIENT  AND  COPYRIGHT  INFRINGEMENT 
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52.227- 10 

52.227- 11 

52.228- 7 

52.230-5 

52.232- 23 

52.233- 1 

52.233- 3 
52.237-3 

52.246- 25 

52.247- 63 

52.249- 5 

52.249- 14 
52.251-1 


FILING  OF  PATENT  APPLICATIONS  - 
CLASSIFIED  SUBJECT  MATTER 

PATENT  RIGHTS  -  RETENTION  BY  THE 
CONTRACTOR  (SHORT  FORM) 

INSURANCE  -  LIABILITY  TO  THIRD 
PERSONS 

COST  ACCOUNTING  STANDARDS  - 
EDUCATIONAL  INSTRUCTIONS 

ALT.  I  -  ASSIGNMENT  OF  CLAIMS 

DISPUTES 

ALT.  I  -  PROTEST  AFTER  AWARD 

CONTINUITY  OF  SERVICES 

LIMITATION  OF  LIABILITY  -  SERVICES 

PREFERENCE  FOR  U.S.  -  FLAG  AIR 
CARRIERS 

TERMINATION  FOR  CONVENIENCE  OF  THE 
GOVERNMENT  (EDUCATIONAL  AND  OTHER 
NONPROFIT  INSTITUTIONS) 

EXCUSABLE  DELAYS 

GOVERNMENT  SUPPLY  SOURCES 
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POD  FAR  CLAUSES 


DESCRIPTION 


252.232-7006  REDUCTIONS  OR  SUSPENSION  OF 

CONTRACT  PAYMENTS  UPON  FINDING  OF 
FRAUD 
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PLANT- WIDE  PREVENTATIVE  MAINTENANCE  AND  MONITORING 


Dr.  Theodore  A.  Bapty 
Dr.  Janos  Sztipanovits 
Jason  Scott 
Brian  Ball 
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Abstract 

This  research  concentrated  on  developing  methods  for  on-line  analysis  monitoring  and 
analysis  of  large-scale  plants.  The  target  application,  monitoring  and  diagnosis  of  a  large  motor 
control  system,  was  chosen  to  focus  the  research  on  a  specific  need  and  to  provide  a  arena  for 
demonstrating  the  technology  through  a  prototype  system.  The  criticality  of  the  motor  system 
to  AEDC,  the  cost  of  an  undetected  fault,  and  the  difficulty  of  diagnosis  guided  the  selection  of 
this  particular  application. 

In  the  course  of  the  work,  several  new  algorithms  were  developed  for  diagnosis  of  the  power 
control  circuits.  The  algorithms  take  advantage  of  Ordered  Binary  Decision  Diagrams  (OBDD) 
to  manage  the  complexity  of  system  diagnosis.  These  algorithms  were  incorporated  into  an 
on-line  monitoring  system  for  connection  directly  to  the  motor  circuit.  The  approach  uses  the 
circuit  diagrams  as  Models.  These  models  are  interpreted  to  generate  the  diagnostic/ monitoring 
system. 
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PLANT- WIDE  PREVENTATIVE  MAINTENANCE  AND  MONITORING 


Dr.  Theodore  A.  Bapty 
Dr.  Janos  Sztipanovits 
Jason  Scott 
Brian  Ball 


INTRODUCTION 

This  research  concentrated  on  developing  methods  for  on-line  analysis  monitoring  and 
analysis  of  large-scale  plants.  The  target  application,  monitoring  and  diagnosis  of  a  large  motor 
control  system,  was  chosen  to  focus  the  research  on  a  specific  need  and  to  provide  a  arena  for 
demonstrating  the  technology  through  a  prototype  system.  The  criticality  of  the  motor  system 
to  AEDC,  the  cost  of  an  undetected  fault,  and  the  difficulty  of  diagnosis  guided  the  selection  of 
this  particular  application. 

The  need  for  this  research  became  apparent  through  participation  in  the  Summer  Research 
Program  (SRP)  at  Arnold  Engineering  Development  Center  (AEDC).  This  report  documents  the 
results  of  a  follow-up  grant  (mini-grant)  awarded  as  an  extension  of  the  SRP. 

In  the  course  of  the  work,  several  new  algorithms  were  developed  for  diagnosis  of  the  power 
control  circuits.  The  algorithms  take  advantage  of  Ordered  Binary  Decision  Diagrams  (OBDD) 
to  manage  the  complexity  of  system  diagnosis.  These  algorithms  were  incorporated  into  an 
on-line  monitoring  system  for  connection  directly  to  the  motor  circuit.  The  approach  uses  the 
circuit  diagrams  as  Models.  These  models  are  interpreted  to  generate  the  diagnostic/monitoring 
system. 


MOTIVATION 

Arnold  Engineering  Development  Center  operates  the  largest,  most  diverse  complex  of 
aerospace  testing  facilities  in  the  world.  These  facilities  range  from  turbine  engine  test  cells, 
to  wind  tunnels,  to  space  chambers,  to  hyperballistic  test  ranges,  to  rocket  test  cells.  Accurate 
simulation  of  conditions  is  a  key  attribute  of  AEDC  testing.  For  turbine  engine  testing,  this 
means  that  air  is  conditioned  to  temperatures,  pressures,  and  mass  flow  rates  encountered  by  a 
full-scale  engine  operating  as  part  of  an  aircraft.  For  wind  tunnels,  the  air  must  be  conditioned 
to  conditions  to  allow  a  scaled-down  model  to  achieve  full  scale  behavior. 

In  order  to  provide  these  simulated  conditions,  a  very  large  complex  of  machinery  is 
required.  For  turbine  engine  testing,  the  Aerospace  System  Test  Facility(ASTF)  uses  6  large 
60,000  horsepower  motors  to  compress  air  through  heaters  and  coolers  to  provide  the  input 
airstream  to  die  test  article.  Another  set  of  12  motors  extracts  the  air  from  the  turbine  engine 
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under  test  to  simulate  altitude  operation.  A  complex  network  of  valves  and  pipes  routes  the  air 
throughout  the  facility  to  provide  a  dynamic  simulation  of  engine  environments. 


Within  AEDC,  there  are  6  major  turbine  engine  test  facilities,  with  multiple  air  supply  and 
exhaust  systems.  The  motors  used  in  these  facilities  range  in  age  from  10  to  50  years.  Many 
of  them  are  irreplaceable.  The  costs  to  repair  and  rewind  one  of  these  motors  is  on  the  order  of 
$3  Million.  More  significantly,  the  downtime  of  the  facility  is  very  expensive,  due  to  the  $20K 
per  hour  value  of  testing  time.  The  impact  of  testing  delays  can  propagate  to  expenses  in  the 
aircraft  system  development  projects. 

One  can  see  the  importance  of  the  correct  functioning  of  the  large  motors  to  the  operations 
at  AEDC.  For  this  reason,  the  motor  control  system  was  chosen  as  a  target  domain. 

The  primary  objective  of  this  effort  is  to  develop  tools  for  the  diagnosis  of  the  large  electric 
motor  systems.  Since  the  motor  itself  is  a  relatively  simple  mechanism,  the  diagnosis  efforts 
focused  on  the  control  system  and  the  environmental  factors,  i.e.  power  supply.  These  factors 
cover  the  majority  of  the  operating  failures  in  the  AEDC  large  motor  complex.  For  problems 
within  the  motor,  indications  are  available  at  the  power  supply.  The  causes  and  indications  of 
failure  of  the  motor  systems  addressed  this  work  included  only  the  electrical  systems.  The 
factors  due  to  mechanical  problems  (bearings,  shaft  imbalance,  system  resonances)  were  not 
considered  within  this  project,  due  to  the  small  scope.  This  is  not  to  say  that  these  factors  are 
not  important:  bearing  wear  can  cause  substantial  problems.  This  information  should  be 
considered  as  a  follow-on  effort. 

The  efforts  will  be  described  in  the  following  sections.  We  begin  with  the  control  system 
diagnosis  efforts,  followed  by  the  power  analysis  results. 


CONTROL  SYSTEM  DIAGNOSTICS 

These  systems  consist  of  a  complex  set  of  interlocks  to  cycle  the  motor  through  its  startup 
phase,  and  to  monitor  for  shutdown  conditions.  The  control  system  determines  when  the  motor 
is  supplied  with  power  and  when  the  power  should  be  removed.  Errors  in  the  control  system 
result  in  improper  motor  functioning.  In  the  best  case,  the  motor  function  is  somewhat  impaired 
or  the  motor  stops.  A  more  serious  malfunction  could  damage  the  motor.  For  these  reasons, 
the  correct  functionality  of  the  control  system  is  critical. 

The  control  system  is  implemented  in  standard  relay  ladder  logic.  In  order  to  manage  the 
complexity  of  this  subsystem,  we  have  chosen  to  use  models  to  represent  its  behavior.  Models 
have  been  used  throughout  history  as  a  way  of  managing  complexity.  The  Measurement  and 
Computing  Systems  Laboratory  at  Vanderbilt  has  been  conducting  research  for  the  past  8  years 
on  using  models  to  describe  complex  systems.  The  Model-Based  approach,  which  we  call  the 
Multigraph  System  Architecture,  has  proven  to  be  effective  in  several  diverse  diagnostics  and 
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monitoring  domains  (see,  for  example,  [1],  [2],  [3]). 

The  essence  of  this  approach  is  as  follows: 

•  Define  a  Modeling  Paradigm.  The  paradigm  defines  the  key  concepts  necessary  to 
capture  information  necessary  to  describe  the  system  to  be  monitored  and  diagnosed. 
The  paradigm  uses  concepts  familiar  to  the  domain,  allowing  the  system  to  be  expressed 
in  a  straightforward,  natural  manner.  Proper  choice  of  paradigm  concepts  is  critical  for 
completeness  of  representation  and  for  ease  of  interpretation. 

•  Implement  a  modeling  editor.  The  editor  is  a  graphical  Computer  Aided  Design  (CAD) 
tool  that  implements  the  modeling  paradigm.  The  model  editor  allows  the  designer  to 
easily  describe  and  input  models  of  the  system. 

•  Implement  an  Interpreter  for  the  models.  The  interpreter  processes  the  models  defined 
in  the  modeling  editor  and  produces  executable  systems.  This  process  is  tailored  to  both 
the  modeling  paradigm  and  the  target  execution  system.  The  interpreter  performs  a 
mapping  from  the  model  paradigm  to  the  execution  system. 

•  Implement  an  Execution  Support  System.  This  software  layer  allows  direct  execution  of 
the  output  of  the  model  interpreter.  The  form  of  the  execution  layer  is  dependent  upon 
the  application  domain.  The  Multigraph  Kernel  provides  a  Large-Grain  Dataflow 
environment  for  implementation  of  parallel  processing  applications.  Diagnosis 
applications  build  a  layer  atop  the  kernel,  to  directly  support  diagnosis  applications.  The 
relay  ladder  logic  of  the  physical  motor  control  system  implements  an  finite  state 
machine.  This  system  is  implemented  in  terms  of  relays  (Coils  and  Contacts),  externally 
driven  contact  closures  and  disconnects,  fuses  and  circuit  breakers,  and  a  network  of 
wires  connecting  these  components. 

Since  the  goal  of  this  project  is  to  diagnose  failures  in  the  control  system,  we  must  consider 
the  types  of  failures  that  occur  within  the  system.  These  are: 

•  Relay  Contact  Closure  Stuck  Open.  Due  to  a  mechanical  failure  of  a  specific  contact 
or  a  failure  of  the  control  coil. 

•  Relay  Contact  Closure  Stuck  Closed.  Due  to  a  mechanical  failure  of  a  specific  contact 
or  a  failure  of  the  control  coil. 

The  most  natural  method  of  modeling  these  systems  is  to  mimic  the  standard  way  of 
representing  the  designs:  A  schematic.  Within  the  modeling  environment  the  following  concepts 
will  be  implemented: 

•  COMPONENTS  drawn  within  a  schematic.  The  supported  components  are: 

D  Externally-Driven  Relay  Contact:  Normally  Open. 

Q  Externally-Driven  Relay  Contact:  Normally  Closed. 

D  Coil-Driven  Relay  Contact:  Normally  Open. 

D  Coil-Driven  Relay  Contact:  Normally  Closed. 
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°  Circuit  Breaker. 

°  Resistor. 

°  Wire. 

All  of  these  components,  with  the  exception  of  the  Wire,  have  two  connection  points. 

Wires  have  only  one  logical  connection  point,  as  the  electrical  potential  is  equivalent 

throughout  a  set  of  connected  wires. 

•  CONNECTIONS  The  components  can  be  connected,  forming  an  electrical  connection. 
Connections  of  wires  form  an  electrical  network.  The  allowable  connections  are: 

°  Wire-to-Wire,  forming  a  net. 

□  Wire-to-Component,  connecting  a  component  to  a  network. 

•  COIL-CONTACT 
RELATIONS  These 
relationships  define  the 
interaction  between 
energizing  a  coil  and 
closure  or  opening  of  a 
contact. 

These  model  concepts  and 
components  are  used  to  draw  the 
control  circuit  within  the 
Model-Builder.  The  information 
necessary  to  build  the  models  is 
available  in  standard  schematic 
diagrams.  These  models  will  be 
interpreted  to  construct  system 
descriptions  suitable  for  the 
diagnostics  algorithms.  These 
algorithms  will  be  described  in  the 
next  section. 


DIAGNOSTICS  ALGORITHMS 

The  goal  of  the  diagnostic 
system  is  to  find  the  faulty 
components  using  observations. 
Since  the  structure  and  expected 
ideal  behavior  of  the  control 
circuit  is  known,  we  choose  a 
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Figure  1:  AFSM  Model  of  the  control  circuit 
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model-based  diagnostic  system  approach  for  solving  the  task.  In  model-based  diagnostics  a  model 
of  the  system  is  used  to  generate  possible  faulty  behaviors.  In  this  process  fault  hypotheses  are 
introduced  in  the  model,  and  the  model-generated  faulty  behavior  is  compared  with  actual 
measurements.  If  the  measured  and  model-generated  behaviors  are  the  same,  the  related  fault 
hypotheses  are  accepted  as  possible  explanations  for  the  fault  symptoms.  Typically,  there  are 
many  fault  hypotheses  that  can  explain  a  particular  fault  symptom.  In  these  cases,  the  diagnostic 
system  selects  the  simplest  possible  explanation,  i.e.  the  one,  which  includes  the  least  number 
of  independent  fault  assumptions. 


1.  System  model 

The  model-based  diagnostic  system  requires  a  system  model,  which  is  able  to  explain  the 
possible  behaviors  in  terms  of  component  characteristics.  The  control  circuit  is  a  relay  logic, 
comprising  switches,  relays  and  contact  points.  This  system  is  modelled  as  an  asynchronous 
finite  state  machine  (ASFM),  as  shown  in  Figure  1.  The  components  of  the  ASFM  model  are 
the  following: 


State  variables:  The  relays  have  two  possible  states,  ’energized’  and  ’not-energized’, 

therefore  their  states  are  modelled  with  Boolean  variables.  Due  to  the 
loops  in  the  switching  network,  the  relay  states  are  considered  to  be  the 
state  variables  of  the  controller  represented  by  the  Boolean  state  vector  j. 
The  state  vector  has  a  ’present  state’,  ^  and  a  ’next  state’  z,  value 
connected  by  an  asynchronous  delay  block.  The  asynchronity  of  the  state 
machine  means  that  only  one  of  the  state  variables  changes  in  the  same 
time. 


Next  State  Logic:  The  relays  are  energized  by  the  potential  difference  appearing  on  the 

’wires’  which  are  connected  to  the  terminals  of  the  relay  coils.  The  ’Next 
State  Logic’  is  represented  by  the  Boolean  expressions  describing  the  next 
state,  of  the  state  variables  in  terms  of  the  states  of  the  ’wires’.  The 
state  of  the  wires  is  characterized  by  the  Boolean  vector  w.  The  state  of 
the  wires  (the  wire  is  ’energized’  or  ’not-energized’)  are  the  directly 
’observable’  variables  in  the  system. 

Switching  Network:  The  state  of  the  wires  is  determined  by  the  state  of  the  relay  contacts,  the 

state  of  the  switches  and  the  topology  of  the  switching  circuit.  The  state 
of  the  switches  and  relay  contacts  can  be  ’open’  or  ’close’.  These  states 
are  modelled  by  the  Boolean  vectors  £,  and 


Switch  Contacts:  The  switch  contacts  are  one  of  the  fault  sources  in  the  control  circuit.  The 

fault  model  of  the  contact  points  is  ’stuck  open’  and  ’stuck  close’.  The 
Switch  Contact  network  models  the  ’real’  switch  contact  points  as  a  logic 
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function  of  the  ideal  switch  contact  state  &  and  a  fault  vector  J.  The  logic 
expression  for  the  fault  model  of  the  ith  contact  point  is  the  following: 

x r(i)  =  xfi)  A  ->X  Jofi )  V  xjc(i), 
x  jofl)  A  xjcfl)  =  0 

where  xr(i)  is  the  ith  real  contact  point,  xfi)  is  the  corresponding  ideal 
contact  point,  x Jofi)  is  the  ’stuck  open’  failure  and  x Jcfi)  is  the  ’stuck 
close’  failure. 

Relay  Contacts:  The  relay  contacts  represent  the  contacts  points  controlled  by  the  state  of 

the  relays.  A  relay  contact  can  be  ’normally  closed’  or  ’normally  open’. 
The  logic  relationship  between  the  s(k)  state  variable  and  the  related, 
c_skfi)  ’normally  closed’  and  c_sk(j)  ’normally  open’  contact  points  is: 

c_sk(j)  =  s(k)  (normally  open) 
c_sk(i)  =  -'sfk)  (normally  closed) 

The  fault  model  of  the  relay  contact  points  is  the  same  as  that  of  the 
switch  contact  points. 

Given  an  initial  state  &  and  an  input  switch  setting  £,  the  state  machine  will  go  through  a 
state  trajectory  and  stabilizes  in  a  new  state  £(&otxJ.  If  the  switch  and  relay  contacts  are  not 
faulty,  the  stabilized  new  state  will  be  the  required  state,  which  e.g.  starts  up  an  electric  motor. 
If  contact  points  are  faulty,  a  the  resulting  new  state  may  be  incorrect,  even  not  stable. 

The  goal  of  the  diagnostic  reasoning  is  to  identify  possible  fault  causes  based  on  the  full, 
or  partial  observation  of  the  state  trajectory. 

2.  Relational  model 

The  main  problem  with  the  diagnostics  of  AFSM  systems  described  above  is  the  very  large 
space  of  possible  behaviors  of  the  system.  Any  attempt  to  perform  diagnostic  reasoning  using 
pattern  matching  seems  to  be  hopeless,  since  in  systems  of  realistic  size  the  fault  patterns  can 
not  even  be  enumerated,  let  alone  used  in  pattern  matching. 
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A  promising  approach  to  solve 
the  state  explosion  problem  is  to 
use  symbolic  reasoning.  In  the 
symbolic  approach,  representations 
focus  not  on  the  individual 
instances  of  variables,  but  on  sets 
and  relations  representing  possible 
states,  fault  hypotheses, 
constraints,  etc.  Sets  and  relations 
can  be  efficiently  represented 
symbolically  using  Ordered  Binary 
Decision  Diagrams  (OBDD)  [1]. 

The  relational  model  of  the 
AFSM  describes  the  system  in 
terms  of  sets  and  relations  as 
shown  in  Figure  2.  The  sets 
represent  all  possible  states  of 
functional  components  of  the 
AFSM.  The  sets  are  constrained 
by  relations  which  are  derived 
from  the  circuit  specification.  For 
example,  the  set  W  includes  all  of 
those  w  E  W  states  of  the  wires 
that  can  be  energized  in  the  same 
time  given  the  operational 
constraints  of  control  circuit.  The 
state  w  is  coded  by  the  Boolean 
vector  w,  i.e.  by  the  states  of  the 
individual  wires. 

The  relations  of  the  AFSM  are  the  following: 

a)  Delay 

A  =  s„fi)  <=>  spfi )  V  i 

b)  State  transition  relation 

R,  =  11,(0)  V  R,(l)  V  ...  V  R,(n-1) 

R,fi)  =  sHfi)  «•  Tt(w)  A  (sji)  ©  spfi))  A  (smQ)  *  s p(j)) 


Figure  2:  Relational  model  of  the  control  circuit 
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where  TJw)  is  a  Boolean  expression  describing  the  next  state  of  the  state  variable  sji)  as 
the  function  of  the  state  of  the  wires,  and  n  is  the  number  of  state  variables  (i.e.  relays). 


c)  Network  relation 

R„  =  RJO)  A  RJI)  A  ...  A  RJm-1) 

R„fi)  =  wfi)  «  Nfa  ,£,) 

where  NJ&,  ,£,)  is  a  Boolean  expression  describing  the  state  of  the  wire  w(i)  in  terms  of  the 
state  of  the  ’real’  switching  contacts  z,  and  the  state  of  the  ’real’  relay  contacts  and  m  is 
the  number  of  wires  in  the  network. 

d)  Switching  contact  relations 

R*  =  RJO)  A  RJI)  A  ...  A  RJp-1) 

RJi)  =  xr(i)  «  (x(i)  A  -ix  Jo  fi)  V  x  Jc(i))  A  (->  xjo(i)  V  ->  xjc(i)) 
where  p  is  the  number  of  switching  contact  points. 

e)  Relay  contact  relations 

Rtf  =  RJO)  A  RJ1)  A  ...  a  RJn-1) 

RJi)  =  RJi.O)  A  RJi.l)  A  ...A  RJi,qrl) 

RJiJ)  =  cjij )*  (c(ij)  A  ~'c Jo(ij)  V  cjc(ij))  a  cjofij)  v  -•  cjcftj)) 

where  cr(ij)  is  the  Boolean  variable  representing  the  jth  ’real’  contact  point  of  the  ith  relay,  and 
qt  is  the  number  of  contact  points  belonging  to  relay  i. 

f)  Relay  constraints 

Re  =  RJO )  A  RJI)  A  ...  A  RJn-1) 

Re(i)  =  RJi.O)  A  RJi.l)  A  ...  A  RcCi,qrl) 

ReCiJ)  =  c(ij)  «  sji) 

where  cfij)  is  the  Boolean  variable  representing  jth  (normally  open)  ’ideal’  contact  point  of 
the  ith  relay. 

3.  Prediction  and  diagnostics 
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Having  the  relational  model  defined,  the  model  can  be  used  for  prediction  and  diagnostics. 

In  prediction,  assuming  a  set  of  initial  states  spESp,  an  input  switch  position  setting  x,  and 
a  given  fault  scenario^  and  fx  the  set  of  next  states  s„E  S„  is  calculated.  The  1-step  prediction 
function  is  defined  as  follows: 

P(SpjcfJJ  =  {sH  |  3  sp,c,cr,w[ (SpGSp)  A  ( (sp,c) ERJ  A  ( (fc,c,cr) E A  ((x/xyxr) E R^)  A 
A  ((xncr,w) ERJ  A  ( w.sj ER, ]} 

Given  the  1-step  predictor  function,  reachability  and  stability  analysis  can  be  performed. 

In  diagnostics,  the  sets  in  the  relational  model  are  propagated  backward.  It  is  important  to 
note  that  the  relational  model  is  not  changing,  only  the  propagation  mechanism.  This  time  we 
describe  the  only  the  batch  diagnostic  method.  In  batch  diagnostics,  the  operation  of  the  control 
system  is  observed,  and  the  results  are  collected  in  the  \m(D,  i n(2),  ...  m(ij ]  measurement 
vector  sequence.  The  components  of  the  m(i)  vector  are  the  states  of  the  observed  wires  at  time 
i.  The  diagnostics  function  is  defined  as  a  recursive  function,  which  propagates  backward  the 
states  and  fault  hypotheses  along  the  fully  or  partially  measured  values  of  the  system  trajectory. 
The  1-step  backward  propagation  function  calculates  the  previous  state  set,  spESp,  the  previous 
wire  state  set  wpE  Wp,  and  the  previous  fault  hypotheses  sets,  f^EF^  and  f^EF^,  given  the 
current  state  set,  snESm,  current  fault  hypotheses  sets,  f^EF^  and  f^EF^,  current  wire  state 
set  wnE  W„  and  measurement  vector  mfi): 

B(SH,F„,F„,Wms,m(i))  =  {sp,  3  sm,c,crfat,fm,wJ(snESJ  A(faEFJ  A(fmEFJ  A 

A  fw.E  WJ  a  (fcpEFJ  A  (fvEFJ  A  ((sp,c)ERJ  A  (( fm,c,cr)E 
ERj)  A  ((x/^JER^  A  ((xr,cr,w)ERJ  A  (( w,sJERt)(m(i)E 
£WJ]} 

Both  the  forward  and  backward  propagation  functions  can  be  calculated  symbolically  using 
the  OBDD  representation  of  the  sets  and  relations.  The  repeated  application  of  the  forward 
propagation  function  generates  the  reachability  set,  the  repeated  application  of  the  backward 
propagation  function  results  in  a  fault  hypothesis  set. 


DIAGNOSTICS  EXECUTION  ENVIRONMENT 

By  interpreting  the  models,  we  can  generate  specifications  for  an  on-line  runtime  system. 

This  system  allows  the  user  to: 

•  Acquire  data  from  the  physical  control  system.  Data  is  logged  and  time-tagged  whenever 
the  state  of  a  sensor  changes.  An  arbitrary  number  of  sensors  is  supported.  The  current 
number  is  24  or  48.  Sample  rate  is  adjustable.  The  hardware  implementation  of  the  data 
acquisition  system  will  be  discussed  later. 
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•  Visualize  the  current  state  of  the  network.  To  facilitate  the  visualization,  a  circuit 
schematic  diagram  is  drawn  on  the  screen.  The  voltage  state  of  each  wire  is  represented 
by  the  color  of  the  wire:  Red  for  high,  black  for  low.  The  component  states  are  shown 
by  changes  in  the  icons.  Activated  coils  change  color  and  closed  contacts  change  visual 
appearance. 

The  user  interface  allows  the  user  to  zoom  in  on  areas  of  the  schematic  and  to  scroll  around. 

•  Step  through  the  historical  data  to  visualize  system  events  at  a  user-controlled  pace.  This 
allows  the  user  to  recreate  events  in  a  time-scale  visible  to  a  human. 

•  Execute  the  diagnostics  starting  from  any  place  in  the  acquired  data  log.  The  diagnostics 
can  be  executed  in  a  forward  simulation  mode  or  in  a  reverse-time  mode. 

•  Visualize  the  results  of  the  diagnostics.  Suspect  components  are  highlighted  in  the 
schematic  display  along  with  a  verbal  description.  The  user  interface  was  implemented 
in  Microsoft  Visual  C+  +  .  The  system  executes  under  Windows  95,  for  maximum 
flexibility  and  ease  of  training  (most  users  already  know  how  to  use  windows 
applications). 


POWER  SYSTEMS  ANALYSIS 

A  parallel  effort  focused  on  the  conditions  of  the  power  supply  feeding  the  motor.  This 
phase  attempts  to  find  current  and  impending  problems  by  monitoring  the  current  and  voltage 
of  the  main  power  feed  to  the  motor. 

For  a  standard,  problem  free  motor  startup  the  design  parameters  of  the  system  specify  a 
range  of  current  and  voltages.  The  first-level  monitoring  system  watches  the  data  online  for 
deviations  outside  the  design  ranges.  Excursions  are  flagged  and  noted  to  a  history  log. 

Each  motor  startup  profile  is  logged  and  parameterized.  The  parameters  for  startup 
characterization  include  rise-time,  slope,  and  damping  parameters  for  the  startup  transient. 
These  parameters  can  be  plotted  versus  time  to  detect  degradations  in  the  motors. 

The  power  consumption  curve  also  corresponds  to  mechanical  motor  load.  By  watching 
instantaneous  Voltage*Current  load  fluctuations  can  be  observed.  Changing  load  profiles  can 
also  indicate  changes  in  motor  behavior. 

The  current  user  interface  does  not  include  the  power  system  analysis.  Live  data  from  motor 
start-up  was  not  available.  Consequently,  the  algorithms  were  not  able  to  be  validated  within 
the  time  frame  of  this  project. 


DATA  ACQUISITION  SYSTEM 

The  data  acquisition  system  used  off-the-shelf  hardware  where  possible  along  with 
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custom-designed  interface  circuitry. 

The  motor  control  circuits  operate  on  125  VDC.  These  voltages  are  too  high  to  feed  directly 
into  standard  computer  input  devices,  which  operate  on  5  VDC.  Direct  contact  of  125V  at  the 
computer  interface  will  damage  the  circuitry.  To  convert  between  these  voltage  levels  and  to 
isolate  and  protect  the  computer,  a  circuit  was  developed. 

The  circuit  consists  of  an  optoisolator  and  a  buffer.  The  125  VDC  is  connected  to  the  input 
to  the  optoisolator  through  a  current-limiting  resistor.  The  output  of  the  optoisolator  drives  a 
TTL  schmitt  trigger  buffer  to  reduce  noise  from  switch  bounce  and  further  protects  the 
computer. 

The  computer  interface  was  implemented  using  QUA-TECH  digital  I/O  boards,  since  they 
were  available  in  the  lab.  Each  of  these  boards  can  acquire  24  digital  inputs.  Any  number  of 
these  boards  can  be  used  in  a  single  computer,  limited  only  by  expansion  slots. 


CONCLUSIONS/RECOMMENDATIONS 

Due  to  the  late  decision  on  the  target  application,  and  to  delays  in  acquiring  information  from 
the  Air  Force  contact,  the  development  of  the  integrated  diagnostics  system  has  not  yet  been 
completed.  The  work  on  integration  of  the  prototype  will  continue  to  bring  the  prototype  to 
completion.  This  prototype  will  be  installed  at  AEDC  to  evaluate  the  effectiveness  of  the 
approach.  We  expect  to  demonstrate  the  prototype  in  late  January,  1996.  The  results  of  the 
prototype  evaluation  will  be  described  in  an  appendix  to  this  report. 

The  diagnostics  methods  developed  for  the  relay  logic  control  circuit  have  been  proven 
effective  in  managing  complexity.  The  original  approach  was  to  describe  all  possible  faults  and 
system  failures  possible  in  the  operation  of  the  system.  Enumerating  all  of  the  failure  modes 
is  a  tedious,  labor  intensive,  error  prone  process.  The  resulting  fault  tree  would  be  very  large, 
with  no  method  for  assuring  complete  system  coverage.  The  OBDD  approach  for  the  execution 
system  support  automatically  ensures  complete  fault  coverage.  This  approach  also  serves  to 
reduce  diagnosis  algorithm  complexity,  contributing  to  system  scalability.  On  the  sample 
application,  the  system  executes  efficiently.  The  scalability  has  been  indicated  by  gradually 
increasing  the  fidelity  of  the  diagnosis,  with  moderate  increases  system  execution  times. 

The  work  with  the  motor  startup  curve  analysis  has  been  entirely  theoretical,  due  to  the  lack 
of  real  data.  For  health  monitoring,  this  data  can  be  invaluable  in  detecting  changes  in  system 
operation.  Analysis  results  look  promising.  Availability  of  high  speed  signal  processing 
hardware  (see  [1])  demonstrates  the  feasibility  of  performing  this  analysis  full-time  at  very  low 
cost.  We  recommend  that  these  approaches  be  further  developed.  The  electrical  aspect  of  the 
motors  is  just  one  part  of  the  operation  of  a  large  system.  Problems  in  the  mechanical 
implementation  of  the  system  can  also  cause  problems,  such  as  bearing  wear,  shaft  imbalance, 
and  gearbox  resonances.  To  detect  and  analyze  these  types  of  problems,  the  mechanical  system 
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must  be  instrumented  with  accelerometers,  vibration  sensors,  and  possibly  audio  sensors. 
Analysis  of  these  signals  can  help  detect  these  problems  and  determine  the  root  causes  of  the 
fault.  Since  a  failure  in  any  aspect  of  the  system  becomes  apparent  in  all  aspects  of  the  system, 
indications  are  visible  in  all  methods  of  monitoring  the  system.  The  most  accurate  system 
diagnosis  will  consider  all  of  these  failure  indicators.  We  recommend  integrating  all  aspects  into 
a  single  system. 
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BLOCK  COPOLYMERS  AT  INORGANIC  SOLID  SURFACES 


John  R.  Dorgan 
Assistant  Professor 

Department  of  Chemical  Engineering  and  Petroleum  Refining 
Colorado  School  of  Mines 


Abstract 


Polymer  molecules  attached  to  a  solid  surface  are  of  extreme  technological  interest.  Thin  polymer  films 
find  widespread  use  in  many  industries  and  improvements  in  thin  films  will  lead  to  superior  technologies. 

In  this  study,  we  report  the  first  experimental  results  on  the  scaling  characteristics  of  brush-forming 
middle-adsorbing  triblocks.  The  triblocks  used  consist  of  relatively  short  polv(ethylene  oxide)  (PEO)  middle 
blocks  and  much  longer  polystyrene  (PS)  end  blocks.  Adsorption  takes  place  onto  a  well  characterized  silicon 
dioxide  surface  from  toluene  and  ellipsometry  is  used  to  determine  the  adsorbed  amount  We  find  that  the  surface 
density,  ct,  for  all  of  the  copolymers  (  both  those  in  the  symmetric  and  asymmetric  regimes  )  scale  according  to  the 
simple  relationships  proposed  in  the  theory  of  Marques  and  Joanny  i.e.  in  the  symmetric  to  moderately  symmetric 


1 

regime,  o  «  — 


1 

regime,  a  «  — - 

P’ 


,  where  NA  is  the  number  of  segments  in  the  adsorbing  PEO  block  and  in  the  highly  asymmetric 
,  where  p  is  the  ratio  of  the  size  of  the  non-adsorbing  block  to  the  size  of  the  adsorbing  block. 
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BLOCK  COPOLYMERS  AT  INORGANIC  SOLID  SURFACES 


John  R.  Dorgan 


Introduction 


Polymer  coils  terminally  attached  to  a  solid  surface  constitute  an  interface  of  particular  interest.  For  example,  they 
are  used  in  the  stabilization  and  flocculation  of  colloidal  particles  such  as  those  used  as  precursors  to  the  sintering 

of  ceramics  Uses  in  biomedical  technology  include  affinity  chromatography  and  the  enhancement  of 
biocompatibility  of  artificial  implants.  These  materials  should  also  prove  useful  for  electrode  modification. 
Composite  materials  consisting  of  adsorbed  polymer  chains  can  be  said  to  represent  the  ultimate  in  polymer  thin 
films;  a  single  layer  of  polymer  coils  attached  to  a  solid  surface.  Such  thin  polymer  films  are  desirable  for  many 
applications  in  electronics  and  optoelectronics  where  film  properties  must  be  closely  controlled 


Objective 

The  objective  of  the  project  was  to  establish  structure-property  relationships  for  polymer  films  formed  via 
the  adsorption  of  block  copolymers  and  to  test  relevant  theories  of  copolymer  adsorption.  The  specific  case  of 
copolymers  to  be  studied  was  the  middle-adsorbing  triblock  copolymers. 


Background 

For  an  isolated  polymer  coil  in  solution,  the  natural  size  of  the  molecule  to  consider  is  the  Floiy  radius, 
Rf-  In  a  good  solvent,  the  Flory  radius  scales  with  the  number  of  repeat  units  raised  to  the  three-fifths  power;  this 
scaling  behavior  represents  a  trade-off  between  the  entropic  penalty  of  stretching  the  coil  and  its  preference  to 

avoid  self-interaction  and  promote  polymer-solvent  contacts  If  a  polymer  molecule  is  attached  to  a  solid  surface 
in  the  presence  of  other  polymer  molecules,  then  the  interchain  spacing  can  become  less  than  the  Flory  radius.  In 
such  a  situation  the  chains  stretch  away  from  the  surface  and  into  the  solvent  -  such  a  structure  is  known  as  a 

3 

polymer  brush  This  stretching  feature  of  polymeric  materials  makes  them  unique  and  contrasts  them  with  low 

4 

molecular  weight  materials  which  face  different  packing  constraints  at  interfaces. 
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With  copolymers,  it  is  possible  to  have  segments  in  the  polymer  chain  which  adsorb  and  segments  which 
do  not.  In  the  discussion  which  follows  the  usual  notation  is  adopted  and  the  (5  parameter  is  defined  as  the  ratio  of 
the  size  of  the  nonadsorbing  (Buoy)  block  to  the  size  of  the  adsorbing  (Anchor)  block 


P  = 


R. 


R. 


aBN° 

a.N° 


■  =  a'(N,/NA)° 


(  1  ) 


In  Equation  1,  aA  and  represent  the  monomer  sizes  within  the  anchor  and  buoy  blocks,  respectively,  while  NA 
and  Nb  represent  the  number  of  monomers  in  each  of  these  blocks.  The  grafting  density  of  polymer  chains  at  the 
surface  also  proves  to  be  an  important  variable  in  the  description  of  polymer  brushes.  The  grafting  density  cr  is 

simply  related  to  the  interchain  spacing  D  through  a  =  (  1  /  D  )^  and  represents  the  number  of  chains  per  unit 
area.  A  schematic  representation  of  a  polymer  brush  in  which  these  various  definitions  are  illustrated  is  shown  in 
Figure  1. 


Scaling  Theory 

For  the  current  work,  the  scaling  theories  of  interest  are  those  involving  the  adsorption  from  non-selective 
solvents.  Marques  and  Joanny  proposed  three  different  regimes  of  symmetry  for  the  case  of  diblocks.  The  results  of 
their  scaling  treatment  are  shown  in  Table  l.5  However,  as  pointed  out  by  Guzonas  et  al,  the  predicted  scaling 
behavior  in  two  of  the  three  regimes  is  similar  and  hence  they  suggested  a  slightly  different  classification  based  on 
the  value  of  (3;  Copolymers  with  P  >  NA3  are  identified  as  highly  asymmetric  (Tail  regime),  those  with 
1  <  P  <  N“  as  moderately  asymmetric  and  those  with  p  <  1  as  symmetric  (Head  regime).6  The  crossover  at  N°5 
is  indicative  of  crossover  from  the  head  regime  (  in  which  the  size  of  the  adsorbing  block  dictates  the  surface 
coverage)  to  the  tail  regime  (  in  which  the  excluded  volume  interactions  of  the  nonadsorbing  blocks  dictates  the 
surface  coverage). 


1-4 


p 

i<P<na0-5 

na°-s<p<na°-7s 

P>na0-75 

3D  regime 

2D  scmidilutc 

2D  dilute 

a  ~ 

Na*1 

P'1 

i 

P'2 

L~ 

n8naw 

Nb^Na2* 

NBi/5NA2/s 

Table  1 :  Scaling  treatment  of  Marques  and  Joanny  for  the  surface  density  and  the  thickness  of  the  brush  layer  for 
diblocks. 


Hence,  the  relative  sizes  of  the  two  blocks  play  an  important  role  in  determining  the  structure  of  the 
surface  layer.  There  are  two  possible  structues  :  either  a  thick,  “fluffy”,  continuous  adsorbed  layer  or  a  thin  layer 
which  may  be  continuous  or  discontinuous. 


Methods 


The  primary  experimental  tools  used  in  this  investigation  were  Low  Angle  Laser  Light  Scattering 
(LALLS),  ellipsometry  (ELLI)  and  Gel  Permeation  Chromatography  (GPC).  Ellipsometry  is  one  of  the  few 
techniques  which  has  the  capability  to  simultaneously  measure  film  thickness  and  grafting  density.  In  order  to 
appreciate  this  point,  ellipsometry  is  now  briefly  reviewed;  a  more  complete  description  of  the  technique  is  readily 
available  in  the  literature.  7 

With  a  linear  polarizer  and  a  half  wave  plate  it  is  possible  to  generate  a  known  state  of  elliptical 
polarization,  when  reflected  from  a  substrate  the  polarization  of  the  light  is  changed.  The  intensity  of  light  passing 
through  the  analyzer  (a  second  motorized  rotating  linear  polarizer)  is  monitored  using  a  photodiode.  From  a 
Fourier  analysis  of  the  intensity  signal,  the  two  ellipsometric  angles  A  and  T,  are  found.  The  fundamental  relation 
of  ellipsometry  relates  refractive  indices  (  n,. )  and  layer  thicknesses  (  d^  )  of  the  reflecting  surface  to  the  measured 
angles  A  and  T. 


e-tan'F^-Ffo.d, 

R ’  (2) 

where  n^.  and  dj.  refer  to  the  indices  of  refraction  and  thicknesses  of  each  layer  present  (denoted  by  the  subscript  k). 
Measurement  of  the  two  independent  quantities,  T  and  A  usually  allows  for  the  solution  of  two  unknowns,  a  layer 
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thickness,  dj  and  a  refractive  index,  n j.  A  separate  determination  of  and  dj  can  be  difficult  and  sometimes 

impossible  due  to  the  small  differences  in  indices  of  refraction  of  the  adsorbate  and  solution.  A  set  of  ellipsometric 
angles,  T  and  A  can  be  produced  by  a  higher  refractive  index,  nj  and  a  lower  thickness,  dj  as  also  a  lower 

refractive  index,  nj  and  a  higher  thickness,  d  j;  but  the  product  n^dj  is  an  invariant  of  the  adopted  layer  model.8 

For  the  study  of  adsorption,  such  a  layer  thickness  and  index  of  refraction  can  be  found  by  assuming  a 
homogeneous  layer.  The  adsorbed  amount  may  then  be  calculated  as: 

A  =  d,cI  =d1(n,-n0)/(dn/dc)o 

Where  nj  represents  the  index  of  refraction  of  the  adsorbed  layer,  c,  is  the  concentration  in  the  layer  and  dj  the 

thickness,  no  represents  the  index  of  refraction  of  the  polymer  solution  and  (dn/dc)0  is  the  change  in  refractive 
index  of  the  solution  with  concentration  of  the  adsorbing  species.  It  has  been  shown  that  the  adsorbed  amount.  A, 
proves  insensitive  towards  which  type  of  concentration  profile  exists  near  the  surface. 

Knowing  the  adsorbed  amount,  the  grafting  density,  a,  is  calculated  from  Equation  4. 

a(chains/ nm!)  =  ■  A^mg/  m  ^  n  (chains/ mol)  x  l(T'*(mJ  /  nm!)  (4) 

M.  (mg  /mol)  "  v 

The  interchain  spacing  may  then  be  calculated  as: 

D  =  (1/ o)°!  (5) 


Further,  the  chain  spacing  at  which  coils  on  the  surface  first  begin  to  touch  is  taken  as: 

D0VHt  =(7tRBJ)O5-(7tRAB:)0!  (6) 

where  Rb  represents  the  radius  of  gyration  of  the  polystyrene  blocks  in  solution  and  R^b  is  the  radius  of  the 
copolymer  in  solution,  determined  by  using  light  scattering. 

The  ellipsometer  was  modified  to  enable  in-situ  studies  of  the  adsorption  studies.  The  details  of  the 
modification  to  accomodate  the  solution  cell  can  be  found  in  previous  works  by  Dorgan  et  al.9 
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Low  angle  laser  light  scattering  (LALLS,  Wyatt  Technologies  DAWN  B)  was  used  to  measure  weight- 
average  molecular  weights  and  the  radii  of  gyration  of  the  copolymer  in  solution.  The  polydispersity  was 
determined  using  gel  permeation  chromatography  (GPC). 


Experiments 

In  this  study,  we  report  on  the  adsorption  behavior  of  middle-attaching  triblocks.  The  materials  employed 
consist  of  relatively  short  polyethylene  oxide)  (PEO)  blocks  capped  by  relatively  long  polystyrene  (PS)  blocks. 
The  PEO  block  preferentially  adsorbs  to  the  surface  whereas  the  PS  block  remains  dangling  in  solution. 
Adsorption  takes  place  from  toluene,  a  good  solvent,  onto  a  well  characterized  silicon  oxide  surface;  this  is  a  case 
of  adsorption  from  a  non-selective  solvent  onto  a  selective  surface. 

The  copolymers  used  in  this  stud)’  have  the  characteristics  shown  in  Table  2.  The  triblocks  were 
synthesized  at  the  Technical  University  of  Istanbul  by  a  procedure  described  earlier.10  It  is  seen  that  the 
composition  of  the  samples  are  such  that  the  different  regions  of  symmetry  are  covered  (i.e. copolymers  with 
P  <  lie  in  the  symmetric  to  moderately  symmetric  regime  and  polymers  with  p  >  N*5  lie  in  the  highly 
asymmetric  regime).  The  silicon  wafers  are  obtained  from  Silicon  Source  Inc.  (  Phoenix)  and  treated  as 
described  in  previous  works.11  The  film  thickness  of  the  oxide  layer  is  independently  determined  prior  to  the 
adsorption  run  using  ellipsometry.  HPLC  grade  toluene  (Aldrich  Chemicals)  was  used  after  filtering  three  times 
through  0.2  micron  Whatman  filters.  The  adsorption  experiments  are  all  conducted  at  solution  concentrations  of 
1.00  ±  0.01  mg/ml.  Low  angle  laser  light  scattering  (Wyatt  Technologies  DAWN  B)  was  used  to  measure  weight- 
average  molecular  weights  and  the  radii  of  gyration  of  the  copolymer  in  solution;  no  evidence  of  micelle  formation 
is  present  at  the  concentrations  used  in  the  experiments.  Figure  2  shows  a  Zimm  plot  for  the  sample  PS(107k)- 
PEO(20k)-PS(107k).  The  polydispersity  was  determined  using  gel  permeation  chromatography  (GPC)  and  seen  to 
be  less  than  3.0  for  all  samples  (  MW/MN  <3.0  ).  A  specially  modified  rotating  analyzer  ellipsometer  (  Gaertner 
Scientific)  is  used  to  measure  the  adsorbed  amounts. 

The  data  analysis  is  performed  as  described  in  previous  section. 


Results  and  Discussions 


Results  from  the  study  are  given  in  Table  2.  It  is  clearly  seen  that  the  ratio  of  the  interchain  spacings  to 
the  overlap  spacings  (  D  /  Dover  )  are  less  than  unity  demonstrating  that  the  triblock  materials  must  be  stretched 
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away  from  the  adsorbing  surface.  To  our  knowledge,  this  is  the  first  report  of  brush  formation  in  copolymers  of  the 
B-A-B  architecture. 


Material 

P 

na0-5 

D  (A) 

D  /  Dover 

PS(254k)-PEO(54k)-PS(254k) 

1.97 

35.03 

284.80 

0.5551 

PS(59k)-PEO(lk)-PS(59k) 

9.16 

4.7 

138.50 

0.3176 

PS(22k)-PEO(  lk)-PS(22k) 

5.07 

4.7 

68.33 

0.1927 

PS(  1 27k)-PEO(  1  k)-PS(  1 27k) 

14.51 

4.7 

205.70 

0.4552 

PS(  107k)-PEO(20k)-PS(  107k) 

2.13 

21.3 

146.90 

0.2929 

PS(305k)-PEO(20k)-PS(305k) 

3.99 

21.3 

248.05 

0.3942 

PS(  168k)-PEO(  10k)-PS(  168k) 

4.26 

15 

195.98 

0.3242 

PS( ’ 1045k)-PEO(  1 0k)-PS(  1045k) 

12.74 

15 

440.11 

0.3492 

Table  2  :  Triblocks  employed  in  this  study  and  the  interchain  spacing  compared  to  overlap  spacing 

Figure  3  shows  that  the  grafting  density  of  the  symmetric  to  moderately  symmetric  triblocks  scale  with 
the  reciprocal  of  the  head  size  (  1/NA  ).  For  highly  asymmetric  triblocks,  the  surface  density  scales  with  the 
reciprocal  of  the  square  of  the  asymmetry  ratio  (  1/p2 )  as  shown  in  Figure  4.  Both  of  these  results  are  consistent 
with  the  predictions  of  the  scaling  theory  of  Marques  and  Joanny  for  diblock  adsorption  The  B-A-B  triblock 
materials  are  thus  seen  to  behave  the  same  as  diblock  materials  as  far  as  the  effect  of  chain  composition  on  the 
structure  of  the  adsorbed  layers  is  concerned. 


Conclusions 


It  is  seen  from  the  (  D  /  Dover)  values  in  Table  2  that  the  B-A-B  triblocks  are  stretched  and  hence  it  can 
be  concluded  that  middle-adsorbing  triblock  materials  can  form  a  brush  structure.  Such  materials  should  therefore 
be  as  effective  as  A-B  diblock  materials  in  steric  stabilization  applications. 

It  can  be  concluded  that  the  middle-adsorbing  triblock  materials  behave  very  similar  to  diblock  materials. 
The  change  in  composition  of  the  B-A-B  copolymers  produces  similar  changes  in  the  structure  of  the  adsorbed 
layer  as  changes  in  the  A-B  diblock  layer.  The  surface  density  for  the  symmetric  to  moderately  symmetric 
triblocks  scale  with  the  reciprocal  of  the  head  size  (1/NA).  For  highly  asymmetric  triblocks,  we  find  that  the 
surface  density  scales  with  the  reciprocal  of  the  square  of  the  asymmetry  ratio  (1/p2)  .  The  B-A-B  architecture 
shares  an  important  feature  with  the  A-B  diblock,  namely  attachment  at  a  single  location  along  the  chain.  This  is 
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in  contrast  to  the  A-B-A  architecture  where  the  chains  may  attach  at  two  different  locations;  this  later  case  may 
produce  distinctly  different  behavior. 
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P  =  a'(NB/NA)°'6 


a  =  (l/D)2 


Figure  1  :  Schematic  representation  of  physical  parameters  of  block  copolymer  adsorption. 
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R(0) 


Figure  2  :  Zimm  plot  for  PS(  107k)-PEO(20k)-PS(107k) 
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STUDIES  OF  SECOND  HARMONIC  GENERATION  IN  GLASS  WAVEGUIDES 

David  Statman 
Associate  Professor 

Department  of  Physics  and  Department  of  Chemistry 
Allegheny  College 


ABSTRACT 


Germanosilicate  glass  waveguides  were  electrically  poled  using  interdigitated  electrodes.  This  poling  resulted 
in  a  quasi-phasematched  second  order  susceptibility.  Second  harmonic  light  was  observed.  Measured  intensities  were 
one  to  two  orders  of  magnitude  smaller  than  expected  based  on  previous  electric  field  induced  second  harmonic  (EFISH) 
generation  studies  with  the  same  interdigitated  electrodes.  The  dynamics  of  the  EFISH  signals  were  studied. 
Exponentially  decaying  signals  suggest  that  through  charge  migration,  via  dielectric  relaxation  and  electronic  drift,  a 
compromising  field  is  established  during  the  poling.  This  results  in  a  weaker  second  order  susceptibility. 
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STUDIES  OF  SECOND  HARMONIC  GENERATION  IN  GLASS  WAVEGUIDES 


David  Statman 


1.  Introduction. 

Second  harmonic  generation  of  laser  light  in  germane-silicate  glass  fibers  and  waveguides  is  a  well  established 
phenomenon..  It  has  been  snggested  that  glasses  which  macroscopically  possess  inversion  symmetry,  the  third  order 

susceptibility  acts  on  an  internal  electric  held.  This  results  in  the  establishment  ofasecond  order  susceptibility  by  which 

second  harmonic  generation  can  be  accomplished..  While  the  mechanism  of  second  harmonic  generation  via  seeding 

is  still  a  matter  of  investigation,  cutrcnt  modeling  is  based  on  a  photorefiactive  effect  where  a  space  charge  held  is 
established  within  the  glass3. 

It  has  also  been  experimentally  confirmed  mat  second  harmonic  light  can  be  genemted  in  glasses  via  the 
ppltcation  of  an  electric  field.  In  fact,  electric  held  induced  second  harmonic  generation  (EFISH)  has  been 
demonstrated  in  glass  fibers-,  and  planar  glass  waveguides.  By  using  interdigitated  electrodes,  Kashyar-demonstrated 
quasi-phasematebed  SHG  in  germanosilicate  fibers,  and  Weitzman,  Hester,  and  Osterberg.  demonstmted  quasi- 

phascmatchedSHGingemumosilicaur  planar  wavepides.  Theirmeasuredsecondhamronic  signal  variedas  the  square 

Of  the  applied  voltage  across  the  electrodes,  as  expected  from  theory. 

Because  electric  fields  are  capable  of  facilitating  second  harmonic  generation  in  glasses,  there  has  also  been 
significant  effort  to  electrically  pole  glasses.  Myers,  Mukherjee,  and  Brueck-.  as  well  as  Naso,  et  aid.  demonstrated 
poling  of  bulk  glass  by  applying  between  I  kV  and  5  kV  across  a  glass  sample  a,  elevated  temperatures.  Each  reported 
second  order  susceptibilities  on  the  order  of  1  pm/V.  Bergot,  et  ala  were  successful  in  poling  glass  fibers  by  fabricating 
them  with  capillary  electrodes  on  opposite  sides  of  the  core.  By  apply  a  dc  electric  field  of  40  kV/cm,  they  were  able 
to  achieve  a  second  order  susceptibility  of  3  fm/V.  Bergot,  et  al.  found  tot  the  second  order  susceptibility  could  be 
enhanced  20  times  with  CW  excitation  by  Argon-laser  light  at  488  run  during  the  poling  process?.  In  order  to  enhance 
the  efficiency  of  second  harmonic  generation  in  optical  fibers,  Kazansky  et  aid  thermally  poled  with  interdigitated 
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electrodes.  The  periodicity  of  the  electrodes  allowed  for  quasi-phasenratching.  By  applyiog  4.3  IcV  across  the 
electrodes,  they  were  only  able  to  obtain  20  pW  of  second  harmonic  for  200  mW  of  IR  pump.  This  amounts  to  a  second 
order  susceptibility  of  about  5  fm/V,  orders  of  magnitude  less  than  expected.  Suggesting  that  a  source  of  electric  field 
spreading  was  electric  breakdown  in  air,  they  were  able  to  improve  the  amount  of  generated  second  hannonic  40  fold 

by  poling  die  fiber  in  a  vacuum10.  Okada.etal.,  were  successful  in  corona  poling  glass  waveguides11.  They  used  a  novel 

technique  for  phasematching  by  fabricating  a  wedge  waveguide  and  finding  the  waveguide  thickness  where  the  TH 

mode  of  the  ftrndamental  was  phasematched  to  the  TM,  mode  of  the  second  hannonic.  With  5  kV  applied  to  a  2.S 

micrometer  thick  waveguide,  they  achieved  a  second  order  susceptibility  of  0.5  pm/V.  Okada,  et  al.  also  found  that 

between  100  and  30OC  poling  temperature  had  no  effect.  Nasu,  et  al  >,  however,  found  tha,  they  could  no.  pole  glass 
at  temperatures  below  100C. 

In  order  to  understand  the  poling  process  better,  Mukherjee,  Myers,  and  Brueck'1  studied  the  dynamics  of  poling 
at  different  temperatures  for  the  case  of  corona  poling  in  bulk  glass.  The,  found  that  the  poling  and  closed  circuit 
depoling  rates  were  chantceriaed  by  a-  1  eV  acdvation  energy.  For  the  open  circuit  depoling  at  elevated  temperature, 

Mukherjee,  e.  al.  measured  an  activation  energy  of  only  0.4  eV.  They  suggested  that  this  was  due  to  surface  charge 
mobility  rather  than  a  bulk  effect. 

In  this  article,  we  investigate  the  effects  of  charge  mobility  on  electric  field  induced  second  hrumonic 

generation  in  waveguides,  and  discuss  the  impact  that  has  on  various  poling  schemes.  By  measuring  the  rime 

dependence  of  EFISH,  we  show  drat  charge  migmdon  results  in  a  compromising  electric  field  which  mitigates  poling 

in  air  using  interdigitated  electrodes.  This  results  in  smaller  than  desired  second  order  suscepdbilides,  jus,  as  described 

by  Kazansky.  The  experiment  is  presented  in  secdon  2.  In  secdon  3,  a  discussion  of  surface  charge  dynamics  and  its 

effect  on  the  internal  electric  field  is  given.  I,  is  shown  drat  such  analysis  agrees  with  the  experimental  data  reasonably 
well. 
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2.  Experiment 

Thermal  Poling  with  Interdigitated  Electrodes. 

Pianargennanosiiicamwaveg™^ 

PO-ea  a,  200*.  Mterfigimmd  eiecncdes  of  the  type  used  by  ^  ^  ^  ^  ^ 

between  the  electrodes  and  the  glass.  The  electrodes  were  spaced  9  5  , 

P  micrometers  apart.  In  order  for  the  poled  field 

to  be  quasi-phasematched  to  the  TM,  modes  of  . 

™.  modes  me  fundamental  wavetengd,  (,.064  micrometem,  and  dm  second 

7"  <532  "*  “  “  ""  «—  —  —  to  me  preplan  aids  of  dm  hmdamenmi.  » 
detennmed  by  measuring  me  propagadon  consmms,  fc  of  me  iiindanmntal  and  second  hannonic  in  ihe 
waveguide  via  prism  coupling1*.  Cam  was  mhen  to  mabe  sure  me  waveguide  and  eiecnodeswereciean.  Shearing 

warn  applied.  (he  ™,  was  momtored  throughout  me  (miing  pmcess  ,0  *  sum  me  msislance  acmss  me 
electrodes  wasneverlessthana  megaohm.  (he  temper  in  me  oven  was  m^ned  a,  200*  for  one  hour  (he 
waveguide  was  allowed  to  cool  overnight  with  the  voltage  applied. 

Ader  poiing.  me  wavegmde  was  probed  for  SHO  using  a  Q-swiicbed  mode  Med  Nd:YAG  (be 
experimental  setup  is  shown  in  Fig.  ,.  (be  wavegui*  was  mouaied  in  a  prism  coup,  er  similar, OMe  used  by  Ubicb 
and  Torge14.  The  output  from  the  laser  was  directed  through  a  halfwave  plateandpoiarizertoinsurethatme  waveguide 

mode  was  TM.  (be  pdsm  coupier  and  wavegui,  were  mounmd  on  a  m,  mb,  ,n  this  way,  the  angie  of  the 

propagaung  beam  could  be  adjusted  with  respect  to  me  grating  vector  of  tbe  poled  field.  The  propagation  length  along 
e  po  ed  portion  of  the  waveguide  was  about  1  cm.  (be  coupled  IRoutpn,  was  monitored  with  a  detected.  Typicaiiy 
the  thtoughpot  of  IR  was  40%.  This  is  consistent  with  a  1  dB  loss  in  each  of  me  coupling  prisms,  and  a  2  dB  loss  along 
die  length  of  me  waveguide.  For 240 mW  average  input  IR  power,  the  IR  power  in  die  waveguide  varied  from  i90  mW 
own  to  120  mW.  The  SHG  signal  was  measured  with  a  photomultiplier  (PM)  whose  output  was  measured  with  a 

boxcar  averager.  The  laser  provided  the  digger  source.  Before  the  PM  were  placed,  filters  to  insure  that  no  IR  signal 
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was  measured.  In  addition,  a  spike  filter  for  wavelength  532  nm  was  introduced.  This  was  to  insure  that  the  signal 
measured  was  SHG  and  not  spurious  fluorescence.  The  spike  filter  typically  reduces  the  intensity  of  532  nm  light  by 

about  50%.  All  signals  were  tested  by  comparing  the  signal  with  and  without  the  filter.  If  the  ratio  was  less  than  50%, 
the  signal  was  rejected  as  not  being  second  harmonic. 

For  the  two  cases  in  which  the  waveguide  was  poled  with  30  V  applied,  second  harmonic  signal  was  detected 
by  the  PM  tube.  On  the  most  sensitive  setting  of  the  boxcar  averager,  second  harmonic  signal  detected  was  well  above 
the  noise  level.  The  signal  to  noise  ratio  exceeded  30.  The  measured  intensities  were  about  50  to  100  picoWatts.  The 
intensities  were  too  low  to  be  seen  visibly.  When  the  waveguide  was  poled  with  1 00  V  applied,  it  was  expected  that  the 
measured  intensity  should  increase  tenfold.  In  thiscase,  however,  the  input  intensity  was  measured  to  beabout400mW, 
giving  waveguide  IR  powers  from  300  mW  down  to  200  mW.  This  resulted  in  IR  degradation  of  the  waveguide  before 
any  intensity  measurement  could  be  made.  However,  before  complete  degradation  of  the  waveguide,  green  light  was 
seen  visibly  in  the  waveguide.  This  suggests  that  the  green  power  was  at  least  a  nanoWatt. 

The  tilt  angle  of  the  poled  electric  field  grating  vector  with  respect  to  the  input  was  varied  to  determine  the  SHG 
efficiency  as  a  function  of  angle.  It  was  found  that  the  acceptance  angle  was  about  3>,  around  the  31-  angle  that  the 
poling  field  grating  vector  made  with  the  fundamental  at  maximum  SHG. 

Maker  fringe  measurements  were  attempted  on  the  poled  waveguides.  No  signal  was  observed.  This  suggests 
that  either  is  less  than  about  5  fm/V,  the  resolution  of  our  instrumentation,  or  the  quasi-phasematching  of  the  poled 

field  (in  and  out  of  the  waveguide)  along  the  length  of  the  waveguide  resulted  in  a  cancelling  out  of  SHG  signal  through 
destructive  interference. 


Dynamics  of  Electric  Field  Induced  Second  Harmonic  (EFISH). 

In  a  second  set  of  experiments,  the  generated  light  was  measured  in  order  to  compare  SHG  signal  from  an 
applied  field  to  that  of  the  poled  waveguide.  The  measured  EFISH  signal  was  significantly  larger  than  the  measured 
SHG  signal  from  the  poled  waveguide.  It  was  also  found,  however,  that  the  EFISH  signal  was  not  robust,  i.e.,  it 
decayed.  This  suggested  a  possible  explanation  for  the  low  SHG  signals  from  the  poled  waveguide.  In  particular,  if  at 
room  temperature  a  compromising  field  is  established  within  the  waveguide  such  that  the  internal  field  is  much  less  than 
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<to““0n  °f  ^  “  —  *>  —  «  —  *om  electric  breakdown  of  *  we  ^ 
argument*?  Hur  ««!•  ,  can  use  no  such 

oUments.  Our  poling  voltages  never  exceeded  imv  ^ 

exceeded  100  V,  as  compared  with  Kazansky  et  al’s  4^no  v  t  „  • 

continually  monitored  the  cutrent  durino  the  nnf  ' 

•  PO  mg  process  and  saw  no  change.  If  there  went  any  electric  breakdown 
in  air,  it  was  certainly  not  across  the  electrodes  I„fcf„ii„  •  'breakdown 

ctrodes.  in  the  following  set  of  experiments,  the  time  dependence  of  tbeEFigH 

signal  was  measured  in  a  second  waveguide. 

ne  experimental  set-np  was  the  same  as  that  shown  in  Fie  1  I„  „„„„  • 

,  ,  Jg-  1-  In  addition,  mterdigitated  electrodes  were 

eZcfi°id  6  !UrfaCe  0f  ^  WaVeS"ide  “  “  ““  ^ “  WeiSmm' et  a”-  “•  ■*>  ■*—  between  dte 
in  prism^  ^  ^cpmpagation  dtrection  was  again  based  on  thep 's  calculated  from  the  coupling  angles 

guide.  The  propagation  length  for  this  waveguide  was  about  2  cm  tv  , 

Th„  ,  cm.  This  resulted  m  losses  of  6  dB  for  the  IR 

^  """  m°re  *■  measured  from  poling. 

Pig  2  shows  the  time  dependence  of  the  SHG  signai  front  the  EFISH  ei  „  r 

FISH  signal  for  five  different  trials  when  the 

r— —  ""**—■ - - 

decays  over  a  period  ofa  few  minutes  After  •  signal 

whenitwasremoved  , nylons,  d ~ -™ei,herwhe„ , he  Md  was  appiied,, nr 

to  decay  to  zero  Raa  •  ?  sig™,  ako  shri„«i  a  decay,  ta  did  no.  appear 

er,  it  appeared  to  come  to  some  steady  state  value.  Again  when  the  field 
eignai  did  go  „P  are,  then  did  decay  to  zero  linv  •  ,  rem°Ve4  ““ SHG 

~^~zzz-=:-:= 

::: 

1,2  and  5  die  IR  pun,  was  always  pKseol  QearIy  ^  ,R 


had  no  significant  impact  on  these  decays. 

Thisdecaydemonstrates  thedynanhcsofacon^nhsin 

by  which  it  was  established. 

In  Hg.  3  the  decaying  SHG  signal  for  when  the  electrodes  were  removed  are  shown  for  three  different  applied 

voltages.  60V,  80V  and  100 V  While  th<»  chnrt  *■  u  u  • 
t  .  '  e*eshortt™^a™^^ 

to  have  the  same  slope  for  all  three  trials. 

Whether  hie  observed  compromising  Held  is  due  to  hie  development  of  a  surface  space  charge  Held,  or  the 
eveloptnent  of  an  internal  space  charge  Heidis  not  dear  from  the  data  presented.  Wha,  is  clear  homfcedata 

re  an  ltfcely  due  to  charge  migration.  Itisnotclearwhatthe  sourceof  theslow decay  is.  InFigs.4and5  the 


b’(l-e  ) 


zrr+  (1  -  A)e  Tz 


(I)  1-  /I(2a>)  A  e-t 

Et  viwr=An 

Where  E  is  die  dc  internal  field  within  die  waveguide.  Tahle  1  gives  the  htpammeters.  The  first  term  of  Eq  lisof 

e  general  form  of  a  decay  due  to  dielectric  relaxation,  electronic  drift  and  diffusion,  whereas  as  the  second  term  is  a 
simple  exponential  decay. 


3.  Discussion. 


In  senem,  the  complex  amplitude  of  the  second  htumonic  A*in  aplanar  waveguide  is  given  hy* 


(2) 


dA2  a 
dz  + 


|A2=i^fWX<2,Aiei^ 
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where  A,  is  to  complex  ampiitude  of  to  todamental,  z  is  to  proton  ditoi„n,«  is  to  attenuafion  coefflden, 
due  .0  linear  losses,  Ms  .he  propagation  constant  ta  to  waveguide,  determined  by  fc  ^  ^  ^ 

accounts  for  the  overlap  integrals  of  the  waveguiding  modes. 

Wilhiu  the  p„,ed  waveguide,  *.  is  deteimined  by  Ihe  interna,  electric  field  tough  fire  tod  order 
susceptibility,  *».  This  poled  field  changes  die  symmehy  of  to  glass  so  to,  1,  no  longer  has  a  macroscopic  inversion 
cento  tomb,  towing  second  hannonic  generton.  Phase  matching  is  accompiished  with  pedothc  poiing  using 
mterdigitated  electrodes.  hi  our  experiment,  thepolingfieidisgiven  approximately  by  EfSintdJJz).  If  to  poling  process 
is  efficient,  X®  ~  X»>E.sin(_apz),  The  second  hannonic  is  phasematched  with  to  fundamental. 

Within  to  nondepleting  pump  approximation,  Eq.  2  can  be  solved  for  to  second  hannonic  field; 

(3)  a  2ko£  „  ... 

A2=i#f(z)*<)E°A?(l-^) 


For  glass  waveguides  to  tod  older  susceptibility  is  *.  -  10-  mW.  and  to  linear  attenuation  ibr  to  sectmd 
harmonic  is  approximately  a = 0.46  cm*.  Foranappiied  fieidof  iOOV,  a  poied  fields  .0  V/pmisexpected,  That 
suggests  tot  X°>  should  be  appioxhnateiy  1  «».  Cieariy  to  is  much  smaller  to  other  reponed  values  of 
approximately  1  pm/v,  where  to  applied  or  poiing  fieid  is  two  orders  of  magnitude  greater  than  that  applied  here  In 

order  to  achieve  higher  second  order  susceptibilities,  iarger  voltages  must  be  applied.  With  our  setup,  however,  no  more 

than  100  V  could  be  applied.  Otherwise,  the  electrodes  would  be  damaged. 

EFiSH  experiments  suggest  to,  even  wi*  to  smai,  vto  of  *»,  the  second  ha^onic  power  in  to  poled 
waveguide  shouid  be  at  leas,  an  order  of  magnitude  greater  ton  measured.  However,  toe  subnaaoWal,  powers  are 
consistent  with  to  experimental  results  of  Bergot,  et.  al. 8  and  Kazansky,  et.  aly  It  would  appear  as  though  the  poling 

process  ts  not  as  efficient  as  one  would  like.  For  our  experiments,  the  electric  field  spreading  due  to  the  breakdown  of 

air  cannot  be  invoked,  since  the  applied  field  is  only  100kV/cm.  Our  FFTSH  pv™  • 

ur  EFISH  experiments  suggest  that  one  possible 

reason  for  this  is  the  establishment  of  a  compromising  charge  field. 


Consider  a  current  density  being  driven  by  the  drift  and  diffusion 


of  charge, 
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(4) 


j  =  }xpcE-DVpc 


where  the  conductivity,  a  is  given  by  the  product  of  the  charge  mobility,  p  and  the  mobile  or  conductive  charge  density 


pc,  o  —  ppc.  From  continuity,  the  time  dependence  of  the  total  charge  density  can  be  determined; 


(5) 


v.j+|=o 


and  hence,  from  Poisson’s  equation  the  time  dependence  of  the  internal  electric  field; 


(b)  ££()V»E  =  p 

Eqs.  4  -  6  can  be  solved  for  an  internal  electric  field  having  the  form  E  =  E^in(APx).  Since  the  total  charge  density 
includes  bound  charges  as  well  as  conductive  charges,  from  Eq.  6  it  can  be  assumed  that  the  conductive  charge  density 
has  the  form  pc  =  p™  +  <t>££o(V»E+ ApEo),  where  <(>  is  a  dimensionless  parameter  corresponding  to  the  amplitude  of  the 
conductive  charge  density  modulated  by  the  electric  field.  p«,  is  the  background  conductive  charge  density.  If  initially, 
at  t  =  0,  the  internal  electric  field  (within  the  waveguide)  is  Ej  with  an  amplitude  E», ,  then  that  field  decays  to  zero  as 


(7) 


with 


E  =  Ei 


e  Ti 


l+bEoi(l-e  T‘  )(1  +cosA{3x) 


and 


Ti 


PPco 

££o 


_l _ 

+  <j>DA{32 


b  = 


PPco 

(JiAjkcc 


t _ 

-  +  DAp 
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The  diffusion  coefficient,  D,  of  a  charge  is  typically  given  by  the  Einstein  relation  D  =  uk.T/e  b  is  then  given  by 


Eq.  7  can  be  compared  to  the  first  term  in  the  empirical  expression  given  by  Eq.  1. 

It  is  clear  fromEq.  7  that  spatially  harmonic  electric  fields  are  generated  during  the  decay  of  the  internal  field. 
The  phase  matched  component  of  the  internal  field,  then  is  given  by  the  first  term  in  Eq.  1  where  b’  is  approximately 
equal  to  bE*.  For  fields  on  the  order  of  lOV/pm  with  A£  =  2jc/A  =  3.3  x  105  m- ,  the  product  bE.  can  be  as  high  as  1.2 
x  103  when  *  is  very  large  (much  greater  than  unity).  Diffusion  and  drift  dominate  the  decay  of  the  field.  For  4  =  0  the 

product  bEoi  equals  0  and  dielectric  relaxation  dominates  the  decay  of  the  field.  In  both  cases,  large  bE,  and  small  bEft 
the  internal  electric  field  decays  to  zero  via  charge  migration. 

For  glass  the  conductivity  o  - 1  x  10-Vohm-n,  Assuming  the  dielectric  relaxation  time  is  much  smaller  than 
the  diffusion  time,  T,  is  approximately  equal  to  the  dielectric  relaxation  time.  Given  that  the  dielectric  relaxation  time 
is  tji  =  EEo/cr,  x,  is  expected  to  be  about  35  s.  It  is  encouraging  that  while  the  spread  of  fit  values  for  x,  is  not  small,  its 
mean  is  about  35  s.  The  expected  and  measured  values  of  X,  are  certainly  within  an  order  of  magnitude  of  each  other. 
The  mobility,  p,  and  the  conductive  charge  density  Pco  cannot  be  calculated  independent  of*.  Nevertheless,  for  a  mean 
value  of  b*  =bEo,  =5  our  measurements  give  p*  =  b/xAP  =  4.3  x  10-mWs,  andp^eeAp/b  =  24  C/m’.  Assuming 
‘t’  t°  be  °n  ^  order  of  unity  (meaning  that  the  compromising  space  charge  field  is  mostly  due  to  conducting  charges  and 
not  bound  charges)  this  gives  a  charge  mobility  of  4.3  x  Iff-  mWs  and  a  conductive  charge  density  of  24  C/itf.  This 
charge  density  corresponds  to  a  carrier  density  of  1.5  x  KWitf.  These  values  are  within  an  order  of  magnitude  of  the 
expected  values  for  charge  mobility  and  carrier  density  in  glass.  Our  results,  then,  suggest  that  the  compromising  field 

is  established  via  dielectric  relaxation  and  electronic  drift.  Diffusion,  however,  is  insignificant  since  it  is  a  much  slower 

process  and  the  charge  modulation  is  not  expected  to  be  large. 
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4.  Conclusions 


We  report  the  observation  of  second  harmonic  generation  in  poled  glass  waveguides.  Measured  intensities, 
however,  were  at  least  an  order  of  magnitude  less  than  expected.  From  EFTSH  measurements  it  was  seen  that  when  an 
electric  field  is  applied  to  the  waveguide,  the  internal  field  decays  in  a  way  that  is  multiexponential.  We  suggest  that 
this  decay  is  the  result  of  charge  migration  through  dielectric  relaxation  and  electronic  drift.  Diffusion  is  not  considered 
to  be  significant.  This  charge  migration  establishes  a  compromising  electric  field,  which  opposes,  at  least  partially,  the 
applied  field. 

In  the  process  of  poling  the  waveguide,  if  these  compromising  fields  are  established,  then  the  poling  field  can 
be  orders  of  magnitude  less  than  that  applied.  The  result  is  a  much  lower  second  order  susceptibility,  hence  much  lower 
second  harmonic  conversion  efficiency.  It  appears  as  though  as  long  as  there  are  enough  conductive  charges  available 
to  create  the  compromising  field,  poling  will  be  inefficient.  This  would  suggest  that  poling  fields  must  exceed  pc/APeeo, 
Eo  >  pc/ApeSo.by  an  order  of  magnitude  or  more.  Assuming  pc  to  be  between  10  and  100  C/m3,  the  poling  field,  then, 
must  be  greater  than  1 07  V/m  by  an  order  of  magnitude  or  more.  While  it  cannot  be  determined  from  these  experiments 
whether  the  charge  migration  is  surface  or  bulk,  surface  charge  effects  can  be  mitigated  if  the  poling  is  done  in  a  vacuum, 
as  Kazansky  has  done. 

Appendix  A: 

In  the  experimental  section,  the  decay  of  the  internal  electric  field  was  fit  to  Eq.  1.  This  equation  was  based 
on  the  assumption  that  the  dynamics  of  that  component  of  the  electric  field  due  to  charge  migration  (as  opposed  to  dipole 
relaxation)  are  determined  by  drift  and  diffusion.  Because  the  conductivity  is  dependent  on  the  charge  density,  hence 
on  the  internal  electric  field,  the  drift  term  is  proportional  to  the  electric  field;  oE  =  aE2.  A  simplified  form  of  this  can 
be  given  as; 


(A-l) 


dE 

dt 
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The  first  term  in  Eq.  A-l  is  due  to  drift,  and  the  second  term  due  to  dielectric  relaxation  and  diffusion.  This  can  be 
solved  to  give 


(A-2) 


E  _  e  T« 

Ej  l+aEoTi(l-e~) 


which  is  Eq.  1,  where  b  =  aE0X!.  A  more  rigorous  derivation  of  the  electric  field  dynamics  can  be  found  in  Section  3. 
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Table  I.  Fitting  parameters  to  Eq.  1  (see  text) 


Inal 

A 

IiM 

hi 

IzM 

1. 

0.43 

62 

4.2 

1200 

2. 

0.58 

14 

~0 

233 

3. 

0.46 

22 

1.6 

250 

4. 

0.48 

37 

8.5 

180 

5. 

0.56 

43 

3.7 

210 
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Figure  Captions: 

Fig.  1 .  Experimental  set-up  for  waveguiding  and  second  harmonic  generation. 

Fig.  2.  SHG  versus  time:  Trials  1-5.  The  applied  voltage  was  BO  V.  Trials  3  and  5  show  the  EFISH  signal  when 

the  field  was  applied,  trials  1,2,  and  4  show  the  EFISH  signal  after  the  field  had  been  removed. 

Fig.  3.  SHG  versus  time  after  removal  of  electric  field,  for  V  =  100V,  80V,  and  60V. 

Fig.  4.  Internal  electric  field  (normalized)  versus  time  for  the  application  of  80  V. 

Fig.  5.  Internal  electric  field  (normalized)  versus  time  for  the  removal  of  80  V. 
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SHG  Intensity  (a.u.) 
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